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13.  ABSTRACT  (MinimimjMtNMU 

Two  storms  from  the  9  August,  1991  CaPE  case  were  analyzed  in-depth  focusing  -?n 
multiparameter  radar  signature  evolution  over  60  min.  in  coordination  with  24 
aircraft  penetrations  which  provided  particle  image  and  electric  field  data  together 
with  vertical  air  motion,  cloud  water  and  other  state  parameters.  A  total  of  five  * 
discrete  "cells"  were  identified  in  the  two  storms  and  their  life  cycle  fully 
documented.  Collaboration  with  South  Dakota  School  of  Mines  and  University  of 
Alabama  at  Huntsville  has  resulted  in  a  full  integration  of  aircraft  image  and  field 
mill  data  (from  SDSM&T  T-28  aircraft)  with  vertical  air  motion  from  dual-Doppler  wind 
synthesis  (UAH) .  The  cellular  evolution  starts  with  a  warm  rain  phase  where  updrafts 
and  a  very  low  concentration  of  large  drops  dominate  the  cloud.  As  the  supercooled 
drops  rise  in  the  updraft  they  freeze  and  acquire  a  water-coat  possibly  by  collisions 
with  other  liquid  drops.  The  multi-parameter  radar  signatures  clearly  identify  this 
mixed-phase  zone.  The  cloud  thereafter  gets  electrified  which  may  intensify  to 
produce  lightning  depending  on  cloud  vertical  growth,  and  generation  of  updraft/ 
downdrafts. 
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g.  QmdssJianmag 

Two  storms  from  9  August,  1991  of  CaPE  were  analyzed,  the  first  one  being  a  multi-cellular 
complex  and  the  second  one  being  more  compact  with  early  identified  and  persistent 
updraft/downdraft  structures.  Three  discrete  cells  from  storm  1  and  two  cells  from  storm  2  were 
analyzed  in  detail  from  near  first  echo  to  mature  stage.  The  first  emphasis  was  on  identification 
of  multiparameter  radar  signatures  such  as  positive  differential  reflectivity  (Z^  columns  extending 
into  sub-freezing  temperatures  of  the  cloud,  enhanced  signatures  of  linear  depolarization  ratio 
(LDR)  "capping"  the  positive  columns,  and  regions  of  enhanced  X-band  specific  attenuation 
(A3).  New  reliable  algorithms  for  estimating  LDR  and  A,  were  developed  taking  into  consideration 
antenna  errors  in  regions  of  strong  reflectivity  gradients,  propagation  effects  and  other  measurement 
errors.  Time-height  contours  through  the  S  cell  cores  were  constructed  for  reflectivity,  Z^,  LDR 
and  A3  which  clearly  showed  the  microphysical  evolution.  Aircraft  data  were  also  integrated  within 
these  time-height  profiles  after  careful  navigation  to  identify  the  cells  penetrated  by  the  various 
aircraft  (4  in  number)  at  various  times  and  altitudes.  Vertical  winds  from  dual-Doppler  analyses 
by  Dr.  Kevin  Knupp  of  the  University  of  Alabama  at  Huntsville  have  been  integrated  into  the  major 
cell  of  storm  1,  and  is  being  also  performed  for  storm  2.  A  comprehensive  manuscript  is  being 
prepared  dealing  with  the  microphysical,  kinematic  and  electrical  evolution  of  the  two  storms  of  9 
August.  Substantial  effort  was  devoted  to  careful  analysis  of  aircraft  particle  image  data  from  which 
particle  types  and  size  distributions  were  obtained.  Careful  consideration  of  navigation  errors  and 
aircraft  "skin  paints"  enabled  the  construction  of  vertical  profiles  of  multiparameter  radar  data  along 
the  aircraft  tracks.  The  aircraft  field  mill  data  from  the  NCAR  King  Air,  T-28  and  P3  were 
carefully  analyzed  for  self-consistency  whenever  the  tracks  penetrated  the  same  cells  at  different 
altitudes.  A  consistent  picture  is  beginning  to  emerge  relating  first  cloud  electrification  to 
appearance  of  mixed-phase  signatures  in  the  multiparameter  radar  data  of  LDR  and  A3.  However, 
persistent  updraft/downdraft  regions  appear  necessary  to  produce,  in  addition,  electric  fields  strong 
enough  to  produce  lightning  discharges.  While  the  transition  from  early,  warm  cloud  to  mature 
phase  is  rather  dramatically  apparent  from  examination  of  positive  columns  followed  by  LDR 
"cap"  and  enhanced  attenuation  near  -5°C  level,  the  strengthening  of  the  electric  field  leading  to  first 
lightning  is  typically  delayed  by  a  further  5-10  minutes  to  allow  for  charge  separation  to  occur  which 
appears  to  be  substantially  accelerated  by  up/down  vertical  air  motions.  A  careful  analysis  of  the 
dual-Doppler  synthesized  vertical  and  horizontal  winds  is  being  performed  to  understand  the 
microphysical-kinematic-electrical  evolution  and  linkages  therein. 

Dr.  I.  J.  Caylor  has  been  spearheading  analysis  of  multiparameter  radar  signatures  in  known 
cases  of  natural  and  triggered  lightning.  The  orientation  if  ice  crystals  in  regions  of  strong  electrical 
fields  has  been  vividly  captured  via  the  increase  in  S-band  specific  differential  phase  (KDP)  followed 
by  a  reduction  to  near  zero  after  lightning  discharge  (crystals  dis-oriented).  The  cyclical  nature  of 
Kqp  correlated  with  lightning  discharges  confirms  McCormick  and  Hendry’s  experimental 
observations  of  nearly  two  decades  earlier  using  circular  polarization  at  K-band.  Multiparameter 
radar  signatures  of  rocket  triggered  lightning  has  been  analyzed  using  polarimetric  and  dual- 
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frequency  reflectivity  measurements  at  S  and  X-bands.  These  results  are  being  further  analyzed  and 
the  number  of  cases  expanded  so  that  some  generalizations  can  be  arrived  at  for  both  natural  and 
triggered  lightning. 
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Analysis  of  the  Copolar  Correlation  Coefficient  between  Horizontal 
and  Vertical  Polarizations 

L.  Liu,  V.  N.  Bringi,  V.  Chandrasekar,  E.  A.  Mueller,  and  a.  Mudukutore 

Colorado  State  University.  Fort  Collins.  Colorado 
(Manuscript  received  20  May  1993,  in  final  form  19  November  1993) 

ABSTRACT 

Recent  research  has  suggested  that  the  copolar  correlation  coefficient  termed  p*»(0)  can  be  used  to  identify 
large  hail  and  improve  polarization  estimates  of  rainfall.  The  typical  measured  values  of  p*.(0)  at  S  band  vary 
approximately  between  0.8  and  1.0.  For  applications  to  hail  identification,  the  required  accuracy  should  be 
within  ±0.01,  while  for  rainfall  improvement  a  higher  accuracy  is  necessary,  for  example,  within  ±0.001.  The 
statistics  of  the  estimator  of  p*.(0)  using  the  Gaussian  spectrum  approximation  from  both  an  analytical  approach 
and  using  simulations  are  discussed.  The  standard  deviation  and  bias  in  p*(0)  are  computed  as  a  function  of 
number  of  samples,  Doppler  spectral  width,  and  mean  value  of  p«.( 0).  The  effect  of  finite  signal- to-noise  ratio 
and  phase  noise  are  also  studied  using  simulations.  Several  other  estimators  rf  p*,(0)  are  evaluated.  Time  series 
data  collected  with  the  Colorado  State  University  CHILL  radar  are  analyzed  and  compared  with  the  simulations. 
Antenna  pattern  effects  as  they  affect  the  accuracy  of  p«>(0)  ate  also  discussed. 


1.  Introduction 

The  magnitude  of  the  zero-lag  copolar  correlation 
coefficient  between  horizontal-  and  vertical-polarized 
received  signals,  termed  as  pte(0),  has  been  recently 
investigated  for  application  toward  hail  sizing,  im¬ 
proving  polarization  estimates  of  rainfall,  and  detection 
of  the  melting  level  in  both  convective  and  stratiform 
precipitation  (Balakrishnan  and  Zmic  1990;  Jameson 
1989;  Zmic  et  al.  1993;  Illingworth  and  Caylor  1991 ). 
A  conventional  differential  reflectivity  ZDr  radar, 
which  alternately  transmits  horizontal-  (H)  and  ver¬ 
tical-  (V)  polarized  waves  (pulse  repetition  time  is  T,) 
and  receives  the  copolar  signal  (either  coherently  or 
incoherently)  via  a  single  receiver,  can  be  used  to  es¬ 
timate  Pfo(0).  A  good  discussion  of  the  Various  physical 
factors  that  affect  this  correlation  are  given  in  Bala¬ 
krishnan  and  Zmic  ( 1990).  Sachidananda  and  Zmic 
( 1985 )  were  the  first  to  consider  the  effect  of  ptn(0)  on 
the  accuracy  of  the  ZDr  measurement.  By  assuming  a 
Gaussian  spectrum,  the  estimator  for  p**(0)  can  be 
written  as  (Sachidananda  and  Zmic  1989) 

Phv(°)  =  ^2Ts)]0-25  ’  (1) 

where  p(2T,)  =  pu,(2T})  -  pm(2Ts)  (see  Fig.  1 ).  It  is 
now  well  established  that  the  standard  deviations  in 
Zdr  and  differential  propagation  phase  $Dp  are  very 
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sensitive  to  pa>(0).  Balakrishnan  and  Zmic  ( 1990)  have 
given  curves  of  the  standard  deviation  and  bias  of 
Pm>(0)  versus  the  number  of  H-V  sample  pairs  (M) 
with  Doppler  spectrum  width  as  a  parameter.  However, 
they  did  not  provide  a  derivation  other  than  stating 
the  method  used,  that  is,  perturbation  analysis  follow¬ 
ing  Zmic  ( 1977). 

Illingworth  and  Caylor  ( 1991 )  estimated  pa*(0)  from 
time  series  of  reflectivity  data  (Z**  and  Zm)  and  showed 
distributions  of  pad(0)  from  light  rainfall  (drizzle)  and 
a  “bright  band,”  the  mean  values  being  0.997  and 
0.849,  respectively.  They  also  provided  a  simple  em¬ 
pirical  relationship  derived  from  data  for  the  standard 
deviation  of  p*v(0).  For  example,  for  a  mean  p**(0) 
value  of  0.99  in  rainfall,  the  standard  error  is  typically 
about  0.002.  From  theoretical  considerations,  that  is, 
assuming  a  gamma  drop  size  distribution  N(D ) 
=  N0Dm  exp[-(3.67  +  m)/D0]  and  equilibrium  oblate 
shapes,  they  found  that  a  change  in  m  from  0  to  2,  at 
a  fixed  ZDR  value,  corresponds  approximately  to  a 
change  in  pap(0)  by  0.002.  Thus,  very  high  accuracy 
in  p*»(0)  is  needed  if  it  is  to  be  used  to  estimate  m  and 
thereby  improve  polarization  estimates  of  rainfall  (see 
also  Jameson  1989).  For  hail  detection  and/or  sizing, 
or  for  brightband  detection,  the  simulations  of  Bala¬ 
krishnan  and  Zmic  (1990)  suggested  that  reduced  ac¬ 
curacy  ( ~  ±0.01 )  in  pm,(0)  may  be  sufficient 

In  this  paper,  we  derive  the  bias  and  standard  de¬ 
viation  of  Phv(0)  given  in  ( 1 )  using  the  perturbation 
method  (details  are  provided  in  appendix  A).  These 
results  are  compared  with  a  simulation  procedure  first 
introduced  by  Chandrasekar  et  al.  ( 1986).  Several  dif¬ 
ferent  estimators  of  p*«,(0)  are  considered  and  com- 
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FK5.  I.  Alternate  copolar  complex  time  sequence  (H — horizontal;  V — vertical). 


pared  using  simulations.  The  effect  of  signal-to-noise 
ratio  (SNR)  and  phase  noise  on  the  bias  of  p*»(0)  are 
also  considered.  Time  series  data  collected  with  the 
Colorado  State  University  CHILL  radar  are  analyzed 
and  stratified  according  to  the  Doppler  spectral  shapes. 
The  experimental  standard  deviations  in  ZD *,  $dp, 
and  i>hv(0)  are  compared  with  simulations.  The  paper 
ends  with  some  general  conclusions  about  the  achiev¬ 
able  accuracy  of  pav(0). 

2.  Statistical  aspects 

a.  Perturbation  method 

The  perturbation  method  is  commonly  used  to  de¬ 
rive  the  standard  deviation  of  nonlinear  estimators. 
This  method  has  been  used  previously  for  deriving  the 
standard  deviation  of  Doppler  velocity  estimators 
(Zmic  1977).  We  extend  this  method  to  derive  the 
bias  and  standard  deviation  of  pm,(0)  defined  in  ( 1 ). 


Figure  1  shows  a  schematic  of  the  alternately  received 
copolar  *ime  samples,  Hv  and  F2l+, .  The  p*.(  1 )  and 
p(2)  are  defined  as 

|V  | 

P*u(  1  )  =  aFH  «  (2) 

(  2  |/fc,|2  2  IK^.I2)05 

i-i  i-i 

iz'jml., 

p(2)  =  pm(2)  *  TJCi - ,  (3) 

2  \Hn\1 

i-i 

which  is  slightly  different  from  the  estimators  given  by 
Balakrishnan  and  Zmic  ( 1990). 

The  bias  and  variance  of  p**(0)  can  be  expressed  as 
(see  appendix  A  for  details) 


biaslpUO)]  =  £[pto(0)]  -  p*(0)  =  pa»(0)|o.15625  “  0.25  ,  (4) 


.  7  ,.,fwlfti(l)]  .  var[pA*(2)l  „  ,  cov[p*,(  1), p**(2)]l 

var[p-<o)l  -  '**'  )(  +  025  -iU2T-°s  )• 
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and  under  the  Gaussian  spectral  assumption  pM(m) 
and  Rhh(m)  can  be  written  as 

Paa("»)  =  expl-SlmT^iyX)2], 

Rhh(m)  =  Rhh(0)phh(m)  exp(—j4mxvT,/\). 

Note  that  v  and  av  are  the  mean  and  spectral  width  of 
the  Doppler  spectrum,  X  is  the  wavelength,  M  is  the 
number  of  H-V  pairs,  and  p**(0)  is  the  true  mean. 

The  bias  predicted  by  (4)  is  very  small  for  typical 
values  of  A/ used  in  current  research  radars  (bias  is  less 
than  0.001  when  M  =  64,  X/4 Ts  =  26.75  m  s-1 ,  pa»(0) 
=  0.95,  and  av  <  4  m  s-1 ).  Figure  2  shows  the  predicted 
standard  deviation  of  pm>(0)  using  (5)  versus  (a)  M 
and  (b)  p*,(0).  In  (a)  as  -  3  m  s-1  and  p*.( 0)  is  a 
parameter,  while  in  (b)  M=  128  with  o’,  as  a  parameter. 
Figure  2b  shows  that  the  standard  deviation  decreases 
more  or  less  linearly  with  increasing  pm,(0),  a  result 
deduced  empirically  by  Illingworth  and  Caylor  ( 1991 ) 
from  their  data.  Their  experimentally  deduced  standard 
deviations  in  rainfall  with  ZDR  in  the  range  of  0-3  dB 
are  less  than  0.004,  which  lies  between  the  curves  for 
pa»(0)  =  0.98  and  0.995  in  Fig.  2a.  However,  it  should 
be  added  that  their  pte(0)  estimates  were  based  on  time 
series  of  power  samples  using  an  interpolation  algo¬ 
rithm  that  they  showed  to  be  nearly  identical  to  the 
power  sample  equivalent  of  ( 1 ). 

b.  Simulations 

Chandrasekar  et  al.  ( 1986)  have  developed  a  spectral 
domain  procedure  for  simulating  the  bivariate  time 
series  of  the  horizontally  and  vertically  polarized  signals 
given  the  Doppler  spectrum,  p**( 0),  and  ZDr.  Green 
( 1990)  has  further  extended  it  to  the  quadvariate  case, 
that  is,  fully  polarimetric  case.  The  simulation  tech¬ 
nique  has  been  used  previously  to  study  and  model 
fluctuations  in  ZDR  and  differential  propagation  phase, 
$dp  (Chandrasekar  and  Bringi  1988;  Chandrasekar  et 
al.  1990;  Bringi  et  al.  1990).  We  use  it  here  to  study 
the  fluctuations  in  pm>(0)  and  to  compare  it  with  the 
perturbation  method.  It  can  also  be  used  to  study  other 
estimators  of  pm,(0)  and  to  study  the  effects  of  SNR, 
phase  noise,  and  non-Gaussian  spectral  shapes  on  the 
bias  and  variance  of  p*,( 0). 

Figure  3  shows  the  simulation  results  wnth  the  same 
parameters  as  in  Fig.  2  (see  Chandrasekar  et  al.  1986; 


Green  1990).  Note  that  ZDR  =  0  dB  in  the  time  series 
simulations  since  the  bias  and  standard  deviation  of 
Pa*( 0)  are  independent  of  Zdr.  The  bias  is  again  very 
small  ( <0.001 )  but  negative  (i.e.,  estimated  values  are 
less  than  the  true  mean ) .  Illingworth  and  Caylor  (1991) 
noticed  that  a  slight  negative  bias  of  about  0.002  would 
bring  good  agreement  between  their  pm>(0)  versus  Zdr 
data  and  theory  for  rainfall.  We  find  good  agreement 
between  the  perturbation  method  and  the  simulations 
(cf.  Figs.  2  and  3),  which  gives  confidence  in  applying 
the  simulation  method  to  study  other  estimators  of 
phv( 0).  Simulations  also  show  that  the  estimator  for 
pMr(0)  used  by  Balakrishnan  and  Zmic  ( 1990)  gives 
bias  and  standard  deviation  that  are  extremely  close 
to  those  obtained  using  (2)  and  (3).  Since  (2)  and  (3) 
are  computationally  less  intensive  we  prefer  it  to  that 
used  by  Balakrishnan  and  Zmic  ( 1990). 


1)  Other  estimators  of  pm,(0) 


The  estimator  defined  in  ( 1 )  will  be  termed  the  al¬ 
ternate  sample  estimator  or  p  2.(0) .  For  incoherent  ra¬ 
dars,  power  samples  can  be  used  to  estimate  p*p(0). 
Let  P/,(2i)  =  | H2i  |2  and p„(2i  +  1)  =  |  K^+i  |2  be  the 
alternately  received  power  samples  as  shown  in  the 
schematic  of  Fig.  1 .  The  estimate  of  the  copolar  power 
correlation  coefficient  pPkP.(  0)  can  be  written  as 


Pr>»(  0)  = 


Pph,p,(  1 ) 

l/W*(2)]025  - 


(9) 


Then, 

pUO)  =  [p,*.*(0)los.  (10) 

This  follows  from  the  fact  that  pp-  p 2 ,  where  pp  stands 
for  pPh  or  p^  (Doviak  and  Zmic  1984).  The  estimator 
given  by  ( 10)  will  be  termed  the  power  sample  esti¬ 
mator.  This  was  one  of  the  estimators  used  by  Illing¬ 
worth  and  Caylor  ( 1991 ). 

Several  methods  can  be  used  to  estimate  the  missing 
received  signal  in  Fig.  1  (Chandrasekar  et  al.  1993) 
(see  also  Fig.  4)  and  in  effect  construct  “simultaneous” 
samples.  Illingworth  and  Caylor  ( 1991 )  used  an  inter¬ 
polation  method  based  on  fast  Fourier  transforms 
(FFTs)  of  power  samples;  their  time  series  data  were 
based  on  a  dwell  time  of  210  ms.  Here  we  apply  a 
frequency  domain  procedure  (Oppenheim  et  al.  1983) 
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Fig.  2.  Theoretically  approximated  standard  deviation  (stdv)  of 
fikv(0)  vs  (a)  number  of  sample  pairs  M,  with  a,  =  3  m  s'1  and 
Pjw(0)  as  a  parameter;  (b)  true  mean  p*,(0),  with  M  -  128  and  a, 
as  a  parameter. 


to  the  alternate  complex  samples  by  first  substituting 
zeroes  for  the  missing  samples  in  Fig.  1  and  then  taking 
the  Fourier  transform  of  the  newly  formed  H'/V'  time 
series  samples  followed  by  ideal  low-pass  filtering  and 
subsequent  inverse  Fourier  transform  to  get  the  inter¬ 
polated  H'/V  time  series  samples  (see  Fig.  4a).  The 
frequency  domain  procedure  yields  a  correlation  coef¬ 
ficient  defined  as 

2  M 

(\/2M)  1  H,V? 

^(0>  ’  (U) 

where  p*  and  pv  are  sample  mean  powers. 

The  estimators  pSUO),  p*®(0),  and  p£,(0)  can  be 
calculated  using  data  from  conventional  ZD r  radars. 
It  is  perhaps  obvious  that  simultaneous  measurements 
of  H  and  V  samples  is  to  be  preferred  for  estimating 
PAt>(0).  Sachidananda  and  Zmic  ( 1985)  consider  a  ra¬ 
dar  model  that  involves  simultaneous  sampling  of  the 
horizontal  and  vertical  polarized  received  signals  from 
two  receivers  when  the  transmit  polarization  is  alter¬ 
nately  switched  between  +45°  and  -45°  linear  states. 


They  used  this  radar  model  to  derive  the  variance  of 
ZD R  and  compare  it  with  the  alternate  sampling 
scheme.  Jameson  and  Dave  ( 1988)  used  circular  po¬ 
larization  radar  data  to  estimate  pm>(0),  which  they 
term  pL.  Figure  4b  shows  a  schematic  of  the  received 
complex  samples  for  the  simultaneous  sampling  case. 
For  a  given  dwell  time,  this  method  will  give  2  M  sample 
pairs,  whereas  the  alternate  method  gives  M  sample 
pairs.  The  copolar  correlation  coefficient  based  on  si¬ 
multaneous  samples  is  also  given  by  ( 1 1 )  and  is  termed 
PUO). 

The  four  estimators  p*»(0),  p^fO),  p£,(0),  and 
P*s(0)  are  compared  in  Fig.  5  using  the  simulation 
procedure.  Both  bias  and  standard  deviation  are  com¬ 
pared  versus  M assuming  true  pto( 0)  =  0.98,  Gaussian 
spectrum  width  at  =  3  m  s~' ,  and  Nyquist  velocity  X/ 
4 T,  =  26.75  ms-1.  The  magnitude  of  the  bias  in 
P*c(0)  is  larger  than  the  other  estimators.  The  standard 
deviation  curves  show  that  p*,( 0)  is  the  best  as  ex¬ 
pected,  followed  by  p£,( 0),  pS»(0),  and  pJU 0).  For 
conventional  coherent  ZDr  radars,  p£,( 0)  is  ideal  for 
Gaussian  spectral  shapes  as  far  as  minimizing  both  bias 
and  standard  deviation  are  concerned.  Computation 


Fig  .  3.  Simulated  standard  deviation  ( stdv )  of  £*(  0 )  vs  ( a )  number 
of  sample  pairs  M,  with  »,*3m  s'*  and  p*.(  0)  as  a  parameter;  (b) 
true  mean  p*(0),  with  M  -  128  and  a,  as  a  parameter. 
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Fig.  4.  (a)  Interpolation  of  missing  points  (dashed  lines)  of  alternate  copolar  sample  sequences; 
(b)  simultaneous  copolar  complex  time  sequence  (H — horizontal;  V — vertical). 


over  a  range  of  true  p*„(0)  and  av  show  that  accuracies 
of  ±0.01  can  be  achieved  with  M>  50  using  p*,( 0). 

2)  Effects  of  additive  noise 

Bringi  et  al.  ( 1983)  derived  a  simple  expression  for 
the  effective  copolar  correlation  coefficient  pJo(0)  in 
the  presence  of  low  SNR.  They  found  that,  for  simul¬ 
taneous  sampling, 

P*r(0)  =  P*»(0)||  1  +—  j|l  +— )]  (12) 

where  both  SNR  and  Zdr  are  in  linear  units.  The  above 
expression  was  verified  by  Illingworth  and  Caylor 
( 1991 )  using  stratiform  rainfall  data  over  a  long-range 
interval  (to  get  different  SNR).  It  is  not  difficult  to 
derive  the  expression  for  the  effective  pZ>(0)  in  the 
presence  of  noise  for  alternate  sampling,  which  is  given 
by 

(1  \-l/4/  Zno  \~,/2 

1  +  snr)  (1  +  sit) 

(13) 

Figure  6  shows  the  results  of  simulations  for  the  alter¬ 
nate  estimator  where  the  bias  in  pj„(0)is  plotted  versus 
SNR  on  a  decibel  scale.  The  solid  line  shows  the  results 
of  correction  based  on  ( 1 3 )  with  Zdr  =  0  dB.  The  fixed 


parameters  are  p* „(0)  =  0.98,  n,  =  3  ms'1,  and  M 
=  128.  It  is  seen  that  the  correction  reduces  the  bias 
significantly.  However,  to  keep  the  bias  within  ±0.01 
for  the  parameters  considered  in  Fig.  6,  the  SNR  should 
be  greater  than  or  equal  to  20  dB.  Figure  6b  shows  the 
results  of  the  corresponding  standard  deviation  versus 
SNR  for  the  effective  pj£(0)  and  bias-corrected 
(>%,( 0).  Again,  for  good  accuracy  the  SNR  should  ex¬ 
ceed  20  dB. 

3)  Effects  of  phase  noise 

The  pulse-to-pulse  rms  phase  error  a4  of  the  complex 
signal  limits  the  clutter  suppression  that  is  possible  in 
conventional,  single-polarized  Doppler  radar  systems 
(Keeler  and  Passarelli  1990).  Klystron  based  trans¬ 
mitters  can  achieve  better  than  0.5°  rms  phase  error, 
while  coherent-on-receive  magnetron  systems  typically 
have  5°-10°  of  rms  phase  error.  A  well-designed  mag¬ 
netron  system  can  achieve  lower  phase  errors  around 
2°-3°. 

The  simulation  method  offers  a  convenient  way  to 
handle  phase  noise.  A  random  phase  of  zero  mean 
Gaussian  (0,  al)  was  added  to  the  simulated  phase  of 
the  complex  signal  H2i  and  V2t+i  in  Fig.  1.  Figure  7a 
shows  the  simulation  results  for  the  bias  in  pjU 0)  ver¬ 
sus  a4  with  pat(0)  =  0.99,  M  =  128,  X  =  5.5  cm,  and 
T,  -  0.83  m  s-' .  The  magnitude  of  the  bias  increases 
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Fig.  5.  Simulated  (a)  bias  and  (b)  standard  deviation  (stdv)  of 
Pto(0)  based  on  estimators  pU( 0),  pj»( 0),  p£( 0),  and  pUO)  vs 
number  of  sample  pairs  M,  with  p*.(0)  =  0.98,  a,  =  3  m  s_l,  and 
X/47,=  26.75  ms-'. 


dramatically  for  <r4  »  4°.  If  the  bias  is  to  be  within  0.01 
for  the  given  parameters  then  <r*  «£  7°.  It  is  obvious 
that  the  bias  in  the  power  estimator  will  be  independent 
of  a*  and  is  about  0.001  for  the  same  parameters.  Thus, 
depending  on  M,  phv(0),  and  <rv ,  there  will  be  a  critical 
value  of  a4  beyond  which  the  power  estimator  is  pre¬ 
ferred.  This  initial  for  the  parameters  used  in  Fig. 
7a  is  around  3°.  Figure  7b  shows  the  corresponding 
standard  deviation  of  pm,(0)  versus  <rt  with  a,  as  a 
parameter.  Again,  the  standard  deviation  increases 
quite  rapidly  for  <r*  >  5°.  For  coherent-on-receive 
magnetron  systems  with  a4  >  3°  it  is  better  to  use  the 
power  estimator  of  p*«,(0)  using  ( 10)  or  the  interpo¬ 
lation  method  used  by  Illingworth  and  Caylor  (1991). 
This  will  avoid  serious  biases  in  pa»(0). 

4)  Effects  of  non-Gaussian  spectra 

The  simulation  procedure  can  also  be  applied  to  non- 
Gaussian  spectra,  for  example,  double-peaked  spectra. 
Multiple-peaked  spectra  can  be  generated  by  super¬ 
imposing  different  Gaussian  spectra,  each  having  dif¬ 
ferent  mean  velocity  v,  spectral  width  av,  and/or  the 


mean  power  P.  An  example  of  such  a  simulated  power 
spectrum  is  shown  in  Fig.  8  where  the  Nyquist  velocity 
is  25  m  s~'.  The  shape  about  each  peak  (at  D,  =  -6 
m  s“‘  and  =  5  m  s')  is  Gaussian  with  av<  =  2  m  s~\ 
P  i  =  1  mW,  and  ^  =  3ms'l,l’j*2  mW,  respec¬ 
tively. 

Figure  9a  shows  the  bias  in  pj>,(0)  (solid  and  dashed 
lines)  and  p£,(0)  (discrete  points)  as  a  function  of  M 
for  symmetric  double-peaked  spectra  with  fixed 
=  =  3  m  s-1  and  S'  =  v2  -  v  as  a  parameter.  Because 

the  alternate  estimator  piv(0)  is  derived  for  Gaussian 
spectrum,  it  is  not  unexpected  to  note  rather  large  biases 
in  pa«(0)  for  non-Gaussian  spectra.  Since  true  p*>( 0) 
was  chosen  to  be  0.98,  the  estimated  values  will  exceed 
1  forSs*  5  m  s-1.  We  emphasize  that  v  is  not  the  mean 
velocity  calculated  for  the  entire  spectrum,  rather  it  is 
the  center  of  one  of  the  Gaussian  shapes  shown  in  Fig. 
8.  Note  also  the  tendency  for  the  bias  to  increase  with 
M  for  v  =  6  m  s-1 .  It  is  clear  that  the  bias  is  serious 
and  pao(0)  should  not  be  used  for  non-Gaussian  spec¬ 
tra.  In  contrast,  the  bias  in  p£,( 0)  decreases  with  M. 
Figure  9b  shows  the  standard  error  of  p*»( 0)  for  v  =  0, 


(a) 


(b) 


Fig.  6.  Simulated  (a)  bias  and  (b)  standard  deviation  (stdv)  of 
pX,( 0)  versus  signal-to-noise  ratio  (SNR)  before  and  after  the  noise 
conection  using  ( 12),  with  p*.(0)  *  0.98,  a,  -  3  ms'1,  and  M 
-  128. 
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FIg.  7.  Simulated  (a)  bias  and  (b)  standard  deviation  (stdv)  of 
pJL(O)  vs  rms  phase  error  <r„  with  p*.(0)  =  0.98,  M  “  128,  and  a, 
as  parameter. 


5,  and  6  ms"'.  Also  shown  is  the  standard  error  of 
p£,( 0)  for  v  =  6  m  s~‘  (the  curves  for  0  and  5  m  s_1 
are  very  near  the  v  —  6  m  s-1  curve  but  lower  and 
hence  not  shown).  By  comparing  the  bias  and  standard 
error  in  pL(0)  and  p£,(0)  for  non-Gaussian  spectral 
shapes,  it  is  clear  that  p £,(0)  is  preferred. 

3.  Radar  time  series  analysis 

The  CSU  CHILL  radar  is  an  S-band,  dual  linear- 
polarized  Doppler  radar  that  can  be  operated  in  a  time 
series  mode  that  allows  recording  of  the  complex  I ,  Q 
signals  for  every  pulse  at  40  range  resolution  volumes, 
which  permits  an  analysis  of  signal  fluctuations.  The 
radar  and  signal  processor  are  described  by  Mueller 
and  Staggs  ( 1986)  and  Brunkow  and  Lee  ( 1986 ).  Re¬ 
cently,  the  signal  processor  was  programmed  to  com¬ 
pute  and  display  p*»( 0)  in  real  time. 

On  25  June  1992  time  series  data  were  collected  for 
several  minutes  with  the  antenna  stationary  and  point¬ 
ing  at  a  convective  precipitation  shaft  at  an  elevation 
angle  of  2°  and  in  the  range  interval  60-70  km  (altitude 
approximately  2.3  km  AGL,  near  the  melting  level). 


riu.  B.  Mraumica  spectrum  magimuuc 
with  two  peaks  centered  at  *  -6  m  s'1  and  Dj  =■  5ms'.  the 
shape  about  each  peak  is  Gaussian  with  s*1 ,  A  =  1  mW, 

and  »tf*3m  s'1,  Pt  «  2  mW,  respectively. 


5  m  s'*.  The 


(a) 


M  (number  of  sample  pairs) 
(b) 


M  (number  of  sample  pairs) 

FIG.  9.  Comparison  between  pL(0)  and  pL(0)  for  Gaussian  spec¬ 
trum  (5  *  0)  and  double-peak  spectrum  (»i  •  5j  “  5,  6  m  s'1, 
respectively).  Simulated  (a)  bias  and  (b)  standard  deviation  (stdv) 
of  p*,(0)  vs  number  of  sample  pairs  M,  with  p*,( 0)  »  0.98  and  a, 
-3m  s'1. 
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Figure  10  shows  range  profiles  of  (a)  reflectivity  Z*, 
(b)  Phv(0),  and  (c)  Doppler  spectrum  width  er„.  The 
average  of  20  range  profiles  (6  s  of  data)  is  shown  to¬ 
gether  with  95%  confidence  interval  and  2 a  standard 
error  bars  for  each  variable.  The  main  precipitation 
core  (peak  Z*  ~  59  dBZ)  is  clearly  visible  at  a  range 
of  64.5  km.  TTie  main  feature  in  p**(0)  is  a  region  of 
lowered  values  between  65  and  67  km  and  the  very 
large  2<r  value  (0.16)  at  66.6  km.  In  this  same  region 
the  Doppler  spectral  width  is  significantly  enhanced 
(peak  of  4.5  m  s-1 )  with  2 a  value  of  2.  /5  ms-1  at  66.7 
km.  The  mean  velocity  (not  shown  here)  changes  quite 
sharply  from  -2.5  m  s-1  at  65  km  to  about  0.5  m  s-1 
at  67  km,  indicative  of  radial  divergence  near  the  melt¬ 
ing  level.  The  intention  in  this  paper  is  not  to  provide 
a  microphysical  explanation  for  the  data,  rather  our 
interest  is  in  the  properties  of  the  pjUO)  estimator. 
Three  resolution  volumes  centered  at  62.3,  65.7,  and 
66.7  km  were  chosen  for  analysis,  the  corresponding 
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Fig.  10.  Average  of  20  range  profiles  of  (a)  reflectivity  Z», 
(b)  p*.(0),  and  (c)  Doppler  spectrum  width  a..  Vertical  lines 
centered  on  the  mean  values  denote  9$%  confidence  intervals 
and  2a  values.  Data  were  collected  by  the  CSU  CHILL  radar 
on  25  June  1992. 


mean  Z*,  p*»(0),  and  av  being,  respectively,  52  dBZ, 
0.95,  1.5  m  s-1;  53  dBZ,  0.87,  2.75  m  s-1;  and  45 
dBZ,  0.89,  4.3  m  s-1.  The  Doppler  spectrum  shapes 
at  the  resolution  volume  at  62.3,  65.7,  and  66.7  km 
were  modeled  based  on  the  measured  spectra  and  are 
shown  in  Fig.  1 1.  The  shapes  are  characterized  as  (a) 
narrow,  (b)  double-peak,  and  (c)  wide.  Using  these 
modeled  spectra  and  true  p*„(0)  of  0.95  for  (a),  0.87 
for  (b),  and  0.89  for  (c),  the  standard  deviations  of 
Zdr,  4>dp,  and  p£t>(0)  were  calculated  versus  M  using 
the  simulation  method  with  A/4 T,  -  26.75  m  s-1 .  The 
results  are  shown  in  Fig.  12  together  with  the  experi¬ 
mentally  measured  standard  deviations.  The  agreement 
between  measurements  and  simulations  is  quite  good 
as  far  as  standard  errors  are  concerned.  The  effect  of 
Doppler  spectral  shape,  while  not  affecting  standard 
errors  in  ZDr  and  4>Dp,  is  quite  pronounced  in  Fig.  1 2c, 
which  shows  the  standard  error  in  ptv(0).  Figure  13 
compares  the  standard  errors  in  the  three  estimators 
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estimator  for  wide  spectral  shapes.  Finally,  Fig.  14 
shows  the  mean  p*»( 0)  versus  M  for  the  three  esti¬ 
mators  based  on  measured  time  series  data  for  (a)  nar¬ 
row,  (b)  double-peak,  and  (c)  wide  spectral  shapes. 


N  0.6 


data, spectrum  a 
aim. , spectrum  a 
data, spectrum  b 
elm. , spectrum  b 
data, spectrum  c 
sim. .  spectrum  c 
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FK>.  1 1.  Modeled  spectra  based  on  radar  data  time  series  mea¬ 
surements  taken  on  25  June  1992:  (a)  narrow,  (b)  double  peak,  and 
(c)  wide. 


M  (number  of  sample  pairs) 


data, spectrum  a 
sim. , spectrum  a 
data, spectrum  b 
sim.  .spectrum  b 
data, spectrum  c 
sim. , spectrum  c 


M  (number  of  sample  pairs) 


data, spectrum  a 
sim. ,  spectrum  a 
data, spectrum  b 
sim. , spectrum  b 
data, spectrum  c 
sim. , spectrum  c 


p% »(0),  pjic(0),  and  p£>(0)  as  a  function  of  M  using 
the  measured  time  series  data  for  (a)  narrow,  (b)  dou¬ 
ble-peak,  and  (c)  wide  spectral  shapes.  While  the 
agreement  between  pS,(0)  and  p*»(0)  is  very  good  for 
(a)  and  (b),  case  (c)  shows  that  p£,(0)  is  the  preferred 


M  (number  of  sample  pairs) 

Fic.  12.  Standard  deviation  (stdv)  of  (a)  2d*,  (b)  *»,  and  (c) 
p  1,(0)  vs  number  of  sample  pairs  M  based  on  data,  collected  by  CSU 
CHILL  radar,  and  simulations  corresponding  to  the  spectrum  of  Figs. 
1  la,  1  lb,  and  1  lc  respectively. 
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Fig.  13.  (a),  (b),  and  (c)  Standard  deviation  (stdv)  of 
P*.(0)  vs  number  of  sample  pairs  M  based  on  CSU  CHILL 
radar  data  corresponding  to  the  spectrum  of  Figs.  I  la,  I  lb, 
and  1  lc,  respectively,  for  pL(0),  pL(0),  and  pf,( 0). 
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Fig.  14.  (a),  (b),  and  (c)  Mean  of  p*( 0)  vs  number  of 
sample  pairs  M  based  on  CSU  CHILL  radar  data  correspond¬ 
ing  to  the  spectrum  of  Figs.  1  la,  I  lb,  and  1  lc,  respectively, 
for  pUO),  p*.(0),  and  pjj,(0). 
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While  good  convergence  is  obtained  experimentally  for 
cases  (a)  and  (b),  case  (c)  shows  that  mean  pjUO) 
does  not  converge  as  a  function  of  M  and  can  exceed 
unity;  thus,  again  p&,( 0)  is  the  preferred  estimator  in 
agreement  with  the  simulation  results  for  non-Gaussian 
spectral  shapes  discussed  in  section  2b. 

4.  Antenna  pattern  effects 

The  previous  discussion  of  the  accuracy  of  pto( 0) 
implicitly  assumes  an  ideal  antenna.  However,  real  an¬ 
tennas  are  not  ideal  and,  in  particular,  phase  matching 
( and  to  a  lesser  extent  amplitude  matching )  of  the  hor¬ 
izontal  and  vertical  polarized  patterns  within  the  main 
lobe  is  crucial  for  accurate  pa*(0)  measurements.  The 
antenna  differential  phase  $da  is  defined  as  0*  -  4>»* 
where  4>h  and  <t>v  are  the  phase  patterns.  Normally,  a 
precise  reflector  illuminated  by  a  well-designed,  dual- 
polarized  horn  will  achieve  a  nearly  constant  $da 
within  the  main  lobe.  The  $da  pattern  of  the  CSU 
CHILL  radar  was  measured  and  it  showed  a  rather 
significant  change  over  the  main  lobe  (Xiao  et  ai. 
1993).  The  antenna  thus  sets  an  upper  bound  on  the 
measured  phv(0).  If  pA  is  the  antenna  coefficient,  it  is 
given  (neglecting  second  order  cross- polar  errors)  by 

|  J  f*fveJ2*°*dSl\ 

Pa  ~  i  »  t>  v  1/2  >  ( 14) 

[Jfldnjfldn) 

where  ft,  and  fv  are  the  conventional  copolar  power 
patterns.  A  two-dimensional  integration  (5°  X  5°  sector 
around  boresight)  of  the  measured  patterns  yielded 
pA  =  0.99,  which  is  the  upper  limit  on  the  observed 
Pac(0).  In  the  case  of  the  CSU  CHILL  antenna  we 
believe  that  numerous  teardowns  and  reassemblies  of 
the  reflector  has  caused  a  deterioration  in  the  reflector 
accuracy.  In  contrast,  the  measurements  reported  by 
Illingworth  and  Caylor  (1991)  in  very  light  rainfall 
show  the  upper  limit  of  their  antenna  to  be  0.998.  Both 
the  National  Severe  Storm  Laboratory  Cimarron  and 
National  Center  for  Atmospheric  Research  CP-2  radars 
appear  to  have  the  upper  limit  of  pa,(0)  limited  by  the 
antenna  to  0.985-0.990  (Balakrishnan  and  Zmic  1990; 
Liu  et  al.  1993). 

5.  Conclusions 

This  study  was  motivated  by  the  findings  of  Illing¬ 
worth  and  Caylor  ( 1991 ),  who  found  that  the  standard 
deviation  of  />**( 0)  obtained  from  time  series  data  did 
not  agree  with  (was  much  smaller  than )  the  theoretical 
predictions  of  Balakrishnan  and  Zmic  ( 1990).  We  have 
derived  the  standard  deviation  and  bias  of  Phv(0)  for 
the  alternate  estimator  using  the  perturbation  approach 
described  by  Balakrishnan  and  Zmic  ( 1990).  We  also 
used  the  simulation  procedure  of  Chandrasekar  et  al. 


( 1 986 )  to  compare  with  the  perturbation  approach  and 
found  good  agreement.  Our  results  also  agree  well  with 
the  experimental  results  of  Illingworth  and  Caylor 
( 1991 ).  A  number  of  other  estimators  of  pm>(0)  were 
analyzed  using  simulations,  that  is,  those  based  on 
power  samples,  simultaneous  samples  of  H  and  V  re¬ 
turns,  and  an  FFI  -based  interpolation  method  used  to 
approximate  simultaneous  sampling.  For  Gaussian 
spectrum,  the  alternate  sample-based  estimator  of 
P*»(0)  gives  smaller  standard  deviations  and  bias  com¬ 
pared  with  the  power  and  FFT  estimators.  For  non- 
Gaussian  spectrum,  the  FFT  method  is  significantly 
better  than  the  alternate  and  power  methods  both  of 
which  give  unacceptably  large  bias.  Other  factors  that 
introduce  bias  are  low  SNR  and  high  phase  noise.  Large 
bias  can  be  introduced  by  low  SNR  for  all  estimators 
when  SNR  <  20  dB.  Simulations  show  that  when  the 
rms  phase  noise  exceeds  about  5°  (e.g.,  magnetron- 
based  systems)  the  power  estimator  of  pa»(0)  should 
be  used  to  avoid  rather  large  biases  that  occur  when 
using  the  alternate  method  based  on  complex  samples. 
A  bias  in  pfo(0)  can  also  be  caused  by  the  antenna  if 
the  phase  difference  pattern  between  H  and  V  polar¬ 
izations  is  not  uniform  across  the  main  beam.  Ampli¬ 
tude  matching  of  the  copolar  patterns,  while  important 
for  accurate  ZM  measurements,  is  not  so  critical  for 
Ph>( 0).  This  depends  on  feedhom  imperfections  and/ 
or  reflector  inaccuracy.  Thus  it  is  exceedingly  difficult 
to  achieve  an  accuracy  level  of  ±0.001  in  pto( 0)  because 
of  these  various  factors.  The  best  results  so  far  have 
been  reported  by  Illingworth  and  Caylor  ( 1991 ),  who 
used  a  dwell  time  of  2 10  ms  in  conditions  likely  to  be 
dominated  by  narrow  spectral  widths  and  Gaussian 
shapes.  Even  then  they  inferred  a  residual  bias  of 
-0.002  in  order  to  bring  their  p^fO)  versus  ZDr  data 
in  better  agreement  with  a  gamma  drop  size  distribu¬ 
tion  model  with  varying  m  parameter.  The  reduced 
accuracy  of  ±0.01  which  is  sufficient  for  hail  detection 
is  readily  attainable  by  using  pjUO)  if  the  spectral  shape 
is  Gaussian  and  by  using  p£,(0)  for  shapes  deviating 
substantially  from  Gaussian.  To  achieve  very  high  ac¬ 
curacies  in  pac(0)  under  all  conditions,  simultaneous 
sampling  is  the  only  alternative  at  the  cost  of  increased 
system  complexity. 
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1.  INTRODUCTION 


In  this  paper,  we  derive  the  bias  and  variance  of  p>,»(0) 
using  the  perturbation  method.  These  results  are  then 
compared  with  simulations.  Several  different  estimators 
of  Phv(O)  are  considered.  Time  series  data  collected  with 
the  CSU-CHILL  radar  are  analyzed. 


2.  STATISTICAL  ASPECTS 

2.1.  Perturbation  Method 

By  assuming  a  Gaussian  spectrum,  the  estimator  for 
Phv(O)  can  be  written  as  ( Sachidananda  &  Zmic,  1989) 

/b..<°)  =  PUi)/[/K2)]°-JS  (1) 


where 


Pa.(l) 

<S(2) 


I  "*va‘+i  I 


[£"i  l*a<  P  YZt  i»fc+i  pio  ! 


£"ii  «*P 


Pha(2)  =  _M 


The  method  described  by  Balakrishnan  &  Zmic  (1990)  is 
used  to  derive  the  bias  and  standard  derivation  of  Pfiv(0) 
defined  in  Eq.  (1).  The  expressions  for  bias  and  variance 
of  Ph.(O)  are  long  and  will  be  given  later  in  a  journal 
paper. 

Fig.l  shows  a  schematic  of  the  alternately  received 
copolar  time  samples,  (Hi,)  and  {Vjj+i}.  Fig.2  shows 
the  predicted  standard  deviation  of  p*»(0)  vs  (a)  M,  and, 
(b)  phv(0).  In  (a)  tr¥  =  3ms-1  and  pa«(0)  is  a  parameter, 
while  in  (b)  M  =  128  with  <r„  as  a  parameter.  Fig.2b 
shows  that  the  standard  deviation  decreases  more  or  less 
linearly  with  increasing  phv(O),  a  fact  deduced  empirically 
by  Illingworth  &  Caylor  (1991)  from  their  data. 

2.2.  Simulations 

We  use  the  simulation  method  ( Chandrasekar  et  al., 
1986)  to  study  the  fluctuations  in  ps«(0)  and  compare  it 
with  the  perturbation  method.  It  can  also  be  used  to 
study  other  estimators  of  pa«(0)  and  the  effects  of  SNR, 
phase  noise  and  non-Gaussian  spectral  shapes  on  the  bias 
and  variance  of  ^hv(O). 

Fig.3  shows  the  standard  deviation  of  vs  (a)  M, 

and,  (b)  pa«(0)  which  can  be  compared  with  Fig.2.  The 
related  bias  is  very  small  (typically  <  0.001).  Excellent 
agreement  is  found  between  the  perturbation  method  and 
the  simulations  which  gives  confidence  in  applying  the 
simulation  method  to  study  other  estimators  of  pa*(0). 

2.3.  Other  Estimators  of  pa»(0) 

The  Estimator  defined  in  Eq.  (1)  will  be  termed  the 
alternate  sample  based  estimator  or  p£,(0).  For  incoher¬ 
ent  radars,  power  samples  can  be  used  to  estimate  p*.( 0). 


v, 


Fig.l  Alternate  co-polar  complex  time  sequence. 
(h  -  horizontal;  v  -  vertical) 


(a)  vs  number  of  sample  pairs  M,  with  oy  =  3  m/s; 

(b)  vs  true  mean  Ps*(0),  with  M  =  128. 


(a)  vs  number  of  sample  pairs  M,  with  o¥  =  3  m/s; 

(b)  vs  true  mean  p/,»( 0),  with  M  =  128. 


(•> 


lb) 


Fig. 4  Simulated  (a)  bias[pf,v(0)]  and  (b)  »tdv[p/,v(0)]  based 
on  different  estimators  vs  number  of  sample  pairs  M, 
with  phu(0)  =  .98  and  <r«  —  3 m/s. 


Let  p/,(2i)  and  p,(2i+l)  be  the  alternately  received  power 
samples  as  shown  in  the  schematic  of  Fig.l.  The  estimate 
of  the  copolar  power  correlation  coefficient  pPk  iP,  (0)  can 
be  written  as, 

Ppt.  .Pm  (0)  =  Ppk.p.  (1)/[pph,pi,  (2)]025 

PhA  o)  =  tPPb.P.(O)]05  (2) 


This  follows  from  the  fact  that  pp  =  p1  where  pp  stands 
for  pPk  or  pp,  ( Doviak  &  Zrnic,  1984).  Estimator  defined 
by  Eq.  (2)  will  be  termed  the  powersample  based  estimate 

pL(°)- 

For  the  simultaneous  sampling  case,  p/.v(0)  can  be  esti¬ 
mated  based  on  2M  pairs  of  H/V  samples  for  a  given  dwell 
time  in  which  alternate  method  gives  M  sample  pairs,  i.e., 

2M 

Pi.(0)=2^£w^*/[PhPv]OS  (3) 

where  pa  and  p„  are  estimates  of  mean  sample  powers. 
Eq.(  3)  is  termed  simutaneous  estimator. 

In  order  to  approximate  the  simultaneous  sampling 
case,  we  apply  the  FFT  technique  on  the  alternate  com¬ 
plex  sample  sequences  for  interpolation  (Oppenheim  & 
Willsky,  1983)  to  estimate  ps«(0)  using  the  Eq.  (3), which 
is  termed  FFT  based  estimator  or  p£v(0)  . 

The  four  estimators  for  p*v( 0)  are  compared  in  Fig. 4 
using  simulations.  The  alternate,  simultaneous  and  FFT 
estimators  are  in  excellent  agreement  as  far  as  bias  and 
standard  deviation  are  concerned.  The  power  estimator 
has  a  small  negative  bias  even  for  large  M  which  is  in 
agreement  with  Illingiiiorth  &  Caylor  (1991).  The  stan¬ 
dard  deviation  is  somewhat  higher  for  the  power  estimator 
as  compared  to  the  other  three. 

3.  TIME  SERIES  MEASUREMENT 

On  25  June  1992  time  series  data  were  collected  by 
the  CSU-CHILL  radar  with  the  antenna  stationary  and 
pointing  at  a  convective  precipitation  shaft  located  in  the 
range  interval  60  -  70  kra.  Three  resolution  volumes  cen¬ 
tered  at  62.3,  65.9  and  66.5  km  were  chosen  for  analy¬ 
sis  which  had  widely  different  Doppler  spectra  as  shown 


in  Fig.5.  Fig.6  shows  the  calculated  standard  deviation 
of  (a)  Zdr ,  (b)  <pd P,  and  (c)  pa»(0)  using  the  alternate 
method  vs  M.  The  three  curves  in  each  panel  correspond 
to  the  three  Doppler  spectra  shown  in  Fig.5  whose  shapes 
can  be  characterized  as,  (a)  narrow,  (b)  double-peak  and 
(c)  wide.  The  discrete  points  in  Fig.6  are  obtained  from 
simulations.  The  standard  deviation  of  pj(¥(0)  shown  in 
Fig.6c  is  very  sensitive  to  spectral  shape,  while  it  is  not 
true  for  Zdr  and  4><tP  as  seen  in  Fig.6a,b. 

Fig. 7a  shows  a  comparison  of  the  standard  deviation 
of  pa«(0)  using  the  three  estimators  p£,(0),  p£„(0)  and 
f>hv(0)  for  the  narrow  spectrum  shape  (Fig.5a).  Similar 
curves  are  shown  in  Fig.7b  for  the  double-peak  spectrum 
(Fig. 5b),  in  Fig.7c  for  the  wide  spectrum  (Fig. 5c).  Corre¬ 
sponding  results  for  estimates  of  the  mean  as  a  function 
of  M  are  given  in  Fig.8.  For  the  narrow  and  double-peak 
spectra  shapes,  the  three  estimators  are  comparable;  for 
good  accuracy  the  power  based  estimator  p£¥(0)  is  not 
preferred.  For  the  wide  spectrum  case  (Fig.7c)  the  al¬ 
ternate  and  power  based  estimators  give  very  high  stan¬ 
dard  deviations  whereas  the  FFT  based  estimator  is  an 
order  of  magnitude  better.  The  mean  estimates  using  the 
three  estimators  are  comparable  for  the  narrow  (Fig. 8a) 
and  double-peak  (Fig.8b)  spectral  shapes.  For  the  wide 
spectrum  (Fig.8c)  only  the  FFT  method  gives  a  conver¬ 
gent  result  while  the  alternate  method  converges  to  values 
slightly  greater  than  unity.  Simulation  results  agree  well 
with  radar  measurements. 

4.  CONCLUSIONS 

We  have  re-derived  the  standard  derivation  and  bias  of 
p/,i>(0)  for  the  alternate  estimator  using  the  perturbation 
approach  of  Balakrishnan  &  Zrnic  (1990).  The  simula¬ 
tion  procedure  of  Chandrasekar  et  al.  (1986)  was  used  to 
compare  with  the  perturbation  approach  and  met  good 
agreement.  Our  results  also  agree  well  with  the  experi¬ 
mental  results  of  Illingworth  &  Caylor  (1991),  and  Caylor 
(personal  communication).  A  number  of  other  estimators 
of  pav( 0)  were  analyzed  using  simulations.  The  accuracy 
in  ph„(0)(±0.001)  for  improving  polarization  estimates  of 
rainfall  is  likely  to  be  achieved  only  by  simultaneous  sam¬ 
pling.  For  hail  detection  the  accuracy  needed  is  an  order 
of  magnitude  lower  (±0.01)  which  can  be  achieved  using 
the  alternate  method  if  the  spectrum  is  Gaussian,  other¬ 
wise  the  FFT  based  method  should  be  used.  Good  agree¬ 
ment  was  found  between  simulations  and  the  CSU-CHILL 
radar  measurements  in  convective  precipitation. 
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Abstract 

Observations  of  lightning  echoes  by  the  NCAR  CP-2  multiparameter  radar  in  cen¬ 
tral  Florida  during  the  1991  CaPE  project  are  presented.  These  observations  are 
unique  in  that  they  are  at  dual  wavelength  and  dual  linear  polarization  and  add 
new  polarimetric  information  to  the  radar  characterization  of  the  lightning  echo. 
An  auxiliary  signal  processor  was  used  during  CaPE  to  make  the  first  polarimetric 
S-band  measurements  of  the  differential  phase  shift  and  copolar  correlation  factor 
of  the  lightning  echo.  The  echoes  were  observed  above  the  melting  level  with  a 
fixed  beam  at  ranges  of  35  to  130  km  and  at  various  heights  between  5  and  12 
km.  Using  time  integrated  radar  data,  statistics  for  the  lightning  echo  are  pre¬ 
sented  for  the  S-  and  X-band  reflectivities  (Zio,  Z3),  differential  reflectivity  (ZDr), 
linear  depolarization  ratio  (LDR),  differential  phase  shift  (<fo p)  and  the  zero  lag 
copolar  cross-correlation  (|/»hv(0)|).  The  statistics  confirm  that  lightning  is  more 
readily  detected  at  S-band  than  X-band  and  indicate  that  |phv(0)|  coupled  with 
Zio  is  a  strong  signature  for  lightning  as  shown  in  Figure  1.  In  addition,  there  is 
a  strong  backscatter  differential  phase  shift  of  up  to  15°  associated  with  the  light¬ 
ning  echo.  A  program  was  carried  out  during  CaPE  to  artificially  trigger  lightning 
using  wire  trailing  rockets.  CP-2  observations  are  presented  in  which  polarimetric 
signatures  are  identified  that  distinguish  between  clouds  where  lightning  was  suc¬ 
cessfully  triggered  and  those  in  which  triggering  failed.  In  general  the  LDR  bright 
band  signature  clearly  indicates  triggering  is  less  successful  in  clouds  composed  of 
low  density  ice  crystals  aloft  and  more  successful  in  clouds  where  LDR  indicates 
higher  concentrations  of  graupel.  In  addition,  the  results  from  an  analysis  of  radar 
time-series  of  a  rocket  triggered  lightning  observed  near  ground  level  are  compared 
to  the  statistics  from  natural  lightning  echoes  in  the  upper  levels  of  thunderstorms. 
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Figure  1.  A  range  profile  of  a  lightning  echo:  a)  showing  enhanced  echo  in  S-band 
reflectivity  (solid)  at  85-86  km  but  no  increase  in  X-band  (dashed).  Note  that 
X-band  power  has  been  attenuated  due  to  heavy  precipitation  at  closer  range  (not 
shown),  b)  showing  approximately  ±10°  backscatter  phase  shift  on  the  <fop  profile 
(solid)  and  the  very  large  decrease  in  |/>hv(0)|  (dashed)  at  85-86  km  range. 
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1.  INTRODUCTION 

The  Convection  and  Precipitation/Electri¬ 
fication  (CaPE)  experiment  was  conducted  in 
east  central  Florida  during  the  period  8  July 
through  18  August  1991.  One  of  the  themes 
of  CaPE  was  coordinated  studies  of  the  co¬ 
evolution  of  the  kinematic,  microphysical,  and 
electrical  properties  of  Florida  thunderstorms. 
We  were  fortunate  to  successfully  perform 
such  coordinated  studies  on  several  occasions, 
using  instrumented  aircraft,  radar,  and  surface 
instrumentation  networks.  We  discuss  here 
observations  of  three  storms  to  illustrate  the 
development  of  precipitation  and  electrifica¬ 
tion  in  loosely-organized  multicellular  Florida 
thunderstorms. 


2.  29  JULY  1991  CASE 

The  South  Dakota  School  of  Mines  and 
Technology  (SDSMT)  T-28  was  launched  on 
this  day  in  anticipation  of  a  boundary  outflow 
collision  to  the  southwest  of  the  multi¬ 
parameter  CP-2  radar.  The  storms  developed 
as  forecast  along  the  collision  zone.  The  T-28 
made  multiple  penetrations  of  one  cell  (cell  A, 
X  -  -17.6  km,  Y  «  -18.6  km  from  the  CP-2) 
at  2137  UTC  and  another  cell  (cell  B, 
X  -  -21.8  km,  Y  -  -12.15  km  from  the  CP-2) 
at  approximately  2154  UTC.  Here  we  follow 
the  evolution  of  cell  A. 
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Cell  A  was  first  detected  by  the  radar  at 
2128  UTC.  By  2130,  its  maximum  echo  inten¬ 
sity  had  increased  to  25  dBz.  The  echo  top  at 
this  time  was  at  6  km  MSL  (-9°C).  The  CP-2 
radar  at  this  time  showed  a  "positive  ZDR 
column”  (indicating  a  region  of  liquid  raindrops) 
extending  up  to  4  km  (4°C). 

Figure  la  depicts  the  time-height  cross- 
section  of  the  cloud  evolution  from  2130  for 
a  period  of  1200  seconds  (20  minutes).  The 
2-dB  ZDR  maximum  height  contour  (indicating 
qualitatively  the  maximum  altitude  at  which 
large  liquid  raindrops  can  be  found  in  the 
cloud)  decreased  from  4  km  (  +  4°C)  gradually 
for  750  seconds,  after  which  moderately 
strong  updraft  (as  measured  by  the  T-28) 
carried  it  back  up  to  4.5  km  (near  0°C). 
Following  this,  the  height  fell  rapidly,  indicating 
fallout  of  the  raindrops.  Figure  1b  combines 
the  reflectivity  and  LDR  history  of  the  cell. 
Regions  of  high  LDR  values  are  observed  after 
800  seconds  at  around  3  km  ( +  10°C)  lying 
well  within  the  high  reflectivity  core,  perhaps 
indicating  a  region  of  mixed  precipitation;  i.e., 
raindrops  with  wet  graupel. 

The  upper  panels  in  Fig.  1  depict  the 
periods  of  T-28  penetrations  and  the  range  of 
vertical  electrical  fields  observed  during  these 
periods.  The  T-28  made  its  first  penetration 
through  this  cell  at  2137  UTC  (420  s  in  Fig.  1) 
at  an  altitude  of  5.5  km  (-6°C).  By  this  time, 
the  peak  reflectivity  had  increased  to  50  dBz, 
as  also  seen  in  Fig.  2,  and  the  top  of  the 
positive  ZDR  column  had  fallen  below  4  km 
( +4°C).  During  this  first  penetration,  the  T-28 
encountered  weak  updrafts  and  downdrafts 
along  with  graupel  and  very  small  amounts  of 
cloud  water.  The  electric  field  was  weak  and 
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Figure  2:  (a)  CAPPI  at  Z  *  5.5  km  at  2137. 
(b)  Vertical  cross-section  depicting  the  Zqr 
contours  at  2137  along  the  path  of  T-28  flight. 
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Figure  3:  (a)  CAPPI  at  Z  «  5.5  km  at  2143. 
(b)  Vertical  cross-section  depicting  the  ZqR 
contours  at  2143  along  the  path  of  T-28  flight. 


mainly  negative,  indicating  positive  charge 
above  or  negative  charge  below  5.5  km. 

The  storm  subsequently  rejuvenated.  By 
2143  UTC  (780  s  in  Fig.  1),  it  had  grown 
considerably  in  size;  the  second  T-28  pene¬ 
tration  began  at  this  time.  Figure  3  depicts  a 
CAPPI  section  of  the  cell  at  aircraft  altitude 
and  the  ZDR  column  along  with  a  vertical 
cross-section  of  ZDR  along  the  aircraft  track. 
The  top  of  the  2-dB  ZDR  column  had  climbed 
to  4.5  km  (0°C)  and  the  echo  top  to  around 
9  km  (-25°C).  A  vertical  radar  reflectivity 
cross-section  along  the  track  of  the  T-28 
shows  peak  reflectivity  of  55  dBz  at  about 
4  km  (+4°C)  and  a  rain  shaft  falling  to  near 
the  ground.  The  maximum  updraft  and  down- 
draft  speeds  encountered  by  the  T-28  during 


this  penetration  at  5.5  km  (-6*0  were  around 
5  m/s  and  the  predominant  hydrometeor  type 
encountered  was  graupel.  The  aircraft's  elec¬ 
tric  field  mills  recorded  positive  vertical  electric 
field  for  this  penetration  with  a  maximum  value 
of  -13  kV/m,  consistent  with  negative  charge 
above  and/or  positive  charge  below  the  flight 
level. 

A  CP-2  scan  at  2146  (960  s  in  Fig.  1), 
right  after  the  aircraft  had  exited  the  cloud, 
shows  precipitation  falling  to  the  ground  (i.e., 
the  maximum  echo  intensity  had  reached  close 
to  the  ground).  At  this  time,  the  top  of  the 
2-dB  ZDR  column  had  fallen  well  below  the 
0°C  level,  indicating  fallout  of  the  main 
precipitation  core  associated  with  a  weakened 
updraft. 
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Organized  updrafts  and  downdrafts  ware  found 
by  aircraft  penetrating  in  the  5  to  6  km  altitude 
range  (-5  to  -9#C)  at  the  time  electrification 
was  first  detected,  but  downdrafts  were  not 
usually  present  before  this.  Measurable  verti¬ 
cal  electric  fields  first  detected  at  these  levels 
were  associated  with  the  downdraft  regions 
containing  graupel  rather  than  with  the  updraft 
regions  containing  cloud  water  with  low 
concentrations  of  rain. 

Takahashi  (19781  found  that,  at  tempera¬ 
tures  above  -10°C,  ice  spheres  became  posi¬ 
tively  charged  during  collisions  with  smaller  ice 
crystals.  However,  at  lower  temperatures  the 
charging  was  a  function  of  cloud  liquid  water 
concentration.  At  either  high  or  low  liquid 
water  concentrations,  ice  spheres  acquired 
positive  charge;  they  acquired  negative 
charges  when  moderate  cloud  water  concen¬ 
trations  were  present.  Gaskell  and  Illingworth 
(1980)  found  similar  results,  indicating  that 
precipitation-size  ice  particles  got  positively 
charged  during  collisions  with  small  ice 
crystals  at  temperatures  in  the  range  -5  to 
-10°C  and  negatively  charged  at  -15  to  -20°C. 

Most  aircraft  penetrations  in  the  three 
storms  discussed  here  were  in  the  -4  to  -10°C 
temperature  range.  The  predominant  sign  of 
the  vertical  electric  field  was  positive,  con¬ 
sistent  with  positive  charge  below  and/or 
negative  charge  above  the  aircraft.  Thus  the  in 
situ  vertical  electric  field  measurements  are 
consistent  with  laboratory  results  for  non- 
inductive  ice  particle  collision  charging  mecha¬ 
nisms.  Particularly,  on  the  last  penetration  of 
the  storm  on  29  July,  when  the  aircraft  was 
near  the  radar  top  of  the  storm  at  5.5  km 
(•6°C),  it  is  almost  certain  that  if  the  pre¬ 
cipitation  particles  at  and  below  this  level  were 
carrying  charge,  this  charge  was  positive. 

As  these  clouds  develop  microphysically, 
they  also  are  evolving  dynamically.  Initially  the 
visible  cloud  is  almost  completely  updraft  with 
compensating  downdrafts  in  the  clear  air 
outside  the  cloud,  so  weak  as  to  be  almost 
undetectable  by  aircraft.  As  rain  forms  by 
coalescence  of  cloud  droplets  within  this 
updraft,  it  is  carried  upward.  It  may  freeze  just 
above  the  freezing  level  (typically  when  the 
convection  is  weak  and/or  disorganized),  or 
perhaps  due  to  ice  ingested  from  a  neighboring 
cloud,  or  it  may  remain  supercooled  to  colder 
temperatures  (typical  in  the  case  of  stronger 
and  more  organized  convection).  As  the  rain 
freezes,  organized  downdrafts  develop  within 


the  cloud,  carrying  the  graupel  grown  upon  the 
frozen-drop  embryos  from  the  5  to  6  km  levels 
to  lower  levels.  The  updraft  also  may  start  to 
collapse  at  this  time.  Another  possibility  is 
that  coalescence  freezing  of  raindrops  could 
result  in  release  of  latent  heat  which  could 
invigorate  the  updraft.  In  the  cases  described 
here,  radar  cloud  top  was  typically  above  8  km 
(-25 °C)  at  the  time  the  2-dB  ZDR  column 
began  to  fall  back  from  the  5  km  (-4°C)  level 
and  electrification  began  to  develop. 

5.  SUMMARY 

All  three  cases  tend  to  show  similar 
signatures.  Following  vertical  development  of 
a  positive  ZDR  column  which  rises  above  the 
freezing  level,  the  raindrops  responsible  for  the 
positive  ZDR  freeze  and  form  a  region  of  high 
LDR  values  aloft.  This  is  followed  by  a  sub¬ 
sequent  increase  in  the  vertical  electric  field 
values  associated  with  downdraft  regions  in 
the  cloud  containing  wet  graupel.  Mixed- 
phase  microphysics  plus  updraft/downdraft 
structures  therefore  appear  to  play  an 
important  role  in  electrification  of  a  cloud. 

A  plausible  explanation  for  the  lack  of 
cloud-to-ground  lightning  for  the  29  July  case 
could  be  that  the  cloud  was  not  very  con¬ 
vective  nor  long-lived.  Its  updrafts  were  weak 
(maximum  measured  updrafts  were  -5  m/s) 
and  short-lived.  The  storms  on  9  August  were 
more  organized  and  grew  to  greater  depths, 
allowing  more  complete  development  of 
precipitation  and  electrification  processes. 
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1.  INTRODUCTION 


2.1  Natural  lightning 


Observations  of  lightning  were  made  in  central  Florida 
during  the  1991  Convective  and  Precipitation/ 
Electrification  experiment  (CaPE)  with  the  NCAR  CP-2 
multiparameter  radar.  Williams  et  al.  (1989)  have 
reported  single  polarization  radar  observations  of 
lightning  at  UHF,  S-band,  and  C-band  while  more 
recently  Krehbiel  et  al.  (1991)  have  employed  a 
wideband  polarimetric  radar  at  X-band.  The  CaPE 
observations  presented  in  this  paper  are  unique  in  that 
they  were  made  simultaneously  at  S-  and  X-band 
wavelengths  with  dual  linear  polarization. 

One  of  the  questions  concerning  the  radar  observation 
of  lightning  deals  with  the  temperature  of  the  plasma 
channel.  Most  recently  Holmes  et  al.  (1980)  have 
interpreted  theii  data  in  terms  of  a  low  temperature 
(underdense)  while  Williams  et  al.  (1989)  have 
supported  the  theory  of  a  hot  (overdense)  plasma.  An 
overdense  plasma  is  a  highly  conducting  target  at 
conventional  radar  wavelengths  and  difficulty  in 
detecting  lightning  at  jhort  wavelengths  is  explained 
by  masking  of  the  lightning  by  the  precipitation  echo. 
The  data  presented  below  address  this  question  and 
the  additional  polarimetric  information  provides 
insight  into  the  spatial  geometry  of  the  lightning  echo. 

2.  OBSERVATIONS 


The  radar  signal  processor  was  configured  to  integrate 
data  for  an  interval  of  128  pulses  (64  horizontal  and  64 
vertical).  The  typical  range  gate  spacing  was  300  m. 

Lightning  echoes  were  observed  in  thunderstorms 
above  the  freezing  level  at  heights  between  5  and  12 
km  and  at  slant  ranges  between  35  and  130  km. 
Figure  1  shows  a  range  profile  with  a  typicai  lightning 
echo  between  555  and  59  km.  The  lightning  is 
particularly  visible  as  a  marked  drop  in  |pHV(0)  |  and 
an  increase  of  Z$  with  a  peak  excursion  of  25  dB  above 
the  precipitation  echo  (58  km).  Note  the  small 
enhancement  of  Zx. 
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The  CP-2  radar  operates  with  a  PRF  of  960  Hz  at  a 
dual  wavelength  of  3.2  and  10.7  centimeters.  The  S- 
band  system  is  coherent  and  provides  estimates  of 
reflectivity  (Z0  and  differential  reflectivity  (ZDR). 
During  CaPE  an  auxiliary  signal  processor  was 
implemented  on  the  S-band  receiver  providing  real¬ 
time  measurements  of  the  differential  phase  shift  («Pop) 
between  the  horizontal  and  vertical  returns  and,  the 
correlation  coefficient  (|Phv(0)|)  between  H  and  V 
returns  (Chandrasekar  et  al.,  1993).  The  X-band 
system  provides  estimates  of  the  reflectivity  (Zx)  and 
linear  depolarization  ratio  (LDR). 

The  data  described  in  this  paper  were  collected  with 
the  CP-2  radar  in  a  manual  mode  with  the  beam 
stationary  in  both  azimuth  and  elevation.  The  X-  and 
S-band  antennas  have  coaxial  beams  with  a  width  of 
about  1®. 


Figure  1.  Range  profile  showing  a  lightning  echo  at 
555  to  59  km.  The  elevation  is  558®  corresponding  to 
height  55  km  above  ground  level. 

Range  profiles  as  in  Figure  1  were  processed  to  isolate 
the  lightning  echo.  The  profiles  were  thresholded  on 
an  S-band  SNR  of  3C  dB  removing  weak,  noisy 
signals.  The  lightning  signature  was  then  isolated  by 
flagging  points  in  the  profile  where  Zj  deviated  by  2.0 
dB  and  |Phv(0)|  deviated  by  0.075  from  their 
respective  means.  The  mean  for  each  parameter  was 
computed  by  averaging  10-15  profiles  (-1.5  s)  prior  to 
the  lightning.  The  flagged  lightning  echoes  were 
extracted  for  further  analysis. 

24  events  spread  over  6  different  days  were  selected 
for  this  paper  from  the  events  logged  during  the  CaPE 
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campaign.  It  is  worth  noting  that  typically  the 
duration  for  each  event  spanned  two  or  more 
integration  cycles  and  extended  over  several 
kilometers  range.  Therefore  155  data  points  were 
available  for  statistic  analysis. 

Using  data  extracted  from  the  lightning  signatures,  the 
distribution  of  (reflectivity  in  units  of  cross-section 
per  unit  volume)  is  shown  in  Figure  2  for  X-  and  S- 
band.  The  corresponding  distributions  Z^,  LDR,  and 
|Phv(0)|  are  shown  in  Figures  3,  4,  and  5. 
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Figure  2.  Reflectivity  distributions  for  S-band  and  X- 
band.  The  mean  for  tj,  is  3x10"*  m3/m3,  while  the 
mean  for  is  3xld7  m2/m). 


Figure  3.  The  differential  reflectivity  distribution  with 
a  mean  value  of  -0.2  dB  and  o=1.2  dB. 


The  differential  phase  shift  upon  backscatter  (8)  was 
estimated  by  subtracting  ipp?  for  each  lightning  echo 
from  the  computed  mean  profile.  A  histogram  of 
the  estimated  6  is  shown  in  Figure  6. 

2.2  Triggered  lightning 

During  the  CaPE  campaign,  NASA  and  French  teams 
engaged  in  triggered  lightning  experiments  using  wire 
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Figure  4.  Distributions  of  the  enhancement  of  the 
copolar  and  cross-polar  X-band  reflectivities  due  to 
lightning. 
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Figure  5.  Distribution  of  the  copolar  correlation  with 
a  mean  of  0.8  and  o  =  0.086. 


Figure  6.  Distribution  of  backscatter  differential  phase 
shift  with  a  mean  of  0.7®  and  a  =  9.5®. 


trailing  rockets.  The  land-based  launch  site  was 
located  approximately  7  km  to  the  south  east  of  CP-2. 
On  25  August  1991,  lightning  was  triggered  at 
175057.130  GMT  at  which  time  the  CP-2  antenna  was 
positioned  over  the  launch  site  at  an  elevation  of  5° 
(height  =  610  m).  The  RP-6  signal  processor  was  in  a 
time-series  collection  mode  with  a  range  gate  spacing 
of  250  m  and  an  integration  interval  of  256  pulses. 

There  was  very  little  precipitation  over  the  launch  site 
at  the  time  of  the  triggering  and  so  precipitation 
masking  is  negligible.  The  radar  echo  had  a  duration 
of  approximately  500  ms,  a  Zp*  of  -4  dB,  Z,  of  32  dBZ. 

3.  DISCUSSION 

If  a  power  law  dependence  on  wavelength  is  assumed 
for  lightning  reflectivity  (t}— X"),  the  value  of  n  can  be 
computed  from  simultaneous  observations  at  two 
wavelengths.  For  each  S-  and  X-band  data  pair  (Fig. 
2)  the  exponent  was  computed.  The  results  in  Figure 
7  show  a  mean  value  near  -2  and  support  the  X'2 
dependence  (bund  by  Williams  et  al.  (1989)  who  used 
UHF,  S-  and  X-band. 

The  time  series  for  the  triggered  lightning  were 
analyzed  in  a  similar  manner  with  the  exponent  being 
computed  for  each  pair  of  dual  wavelength  pulses. 
Surprisingly  the  mean  value  for  n  is  near  +1.5  (Fig.  7). 
A  straightforward  interpretation  of  this  result  is  that 
the  channel  is  underdense.  However  video 
photography  of  the  triggered  lightning  indicates  the 
channel  had  a  large  3  dimensional  curvature  and  so  it 
is  possible  that  the  lightning  was  being  observed  in 
the  antenna  sidelobes.  Further  analysis  of  these  time 
series,  including  a  Doppler  spectral  analysis,  along 
with  ground  based  electrical  data  is  planned. 


Figure  7.  Distribution  of  the  power  law  exponent,  n, 
for  natural  (mean  of  -1.7)  and  triggered  lightning 
(mean  of  15). 


In  general  the  polarimetric  signatures  (Figs.  3-6) 
indicate  the  lightning  target  is  geometrically  complex. 
The  wide  variation  of  Zf*  and  &  with  means  near  zero 
indicate  no  preferred  orientation  of  the  target.  The  2.4 
dBZ  enhancement  of  the  X-band  crosspolar  reflectivity 
is  indicative  of  a  randomly  oriented  target  as  well. 

Analysis  shows  that  decreases  in  |pHV(0)|  are  not 
correlated  with  large  B.  This  could  be  due  to  multiple 
scattering  interactions  between  individual  lightning 
elements  within  the  pulse  volume. 

4.  CONCLUSION 

Multi  para  meter  radar  observations  of  natural  lightning 
during  the  CaPE  project  show  the  following: 

•  The  parameter  |Phv(0)|  decreases  dramatically  for 
the  lightning  echo  and  coupled  with  Z$  is  an  excellent 
polarimetric  signature  for  lightning. 

•  There  is  a  X2  dependence  between  X-band  and  S- 
band  reflectivities  for  lightning  in  precipitation. 

•  X-band  cross-polar  results  indicate  the  lightning 
echo  is  depolarizing.  ZpR  data  show  that  there  is  no 
apparent  channel  orientation.  In  addition  there  is  a 
large  backscatter  differential  phase  shift  of  up  to  ±15° 
which  is  uncorrelated  with  low  |pHV(0)|. 

•  The  only  triggered  lightning  observed  during  CaPE 
displays  an  unusual  signature.  Further  observations 
of  triggered  lightning  would  be  of  great  benefit. 
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1.  INTRODUCTION 

On  9  August,  1991,  convection  began  forming  over  the 
Cape  Canaveral  area  at  16:30  UTC  due  to  interactions  be¬ 
tween  horizontal  convective  rolls  and  river  breezes.  Sur¬ 
face  winds  were  from  the  W/NW  near  the  Cape  area.  The 
15:00  sounding  at  Cape  Canaveral  indicated  light  west¬ 
erly  winds  up  to  6  km  which  then  shifted  to  E/SE  winds 
above  7.5  km.  The  0*C  level  was  at  4.6  km,  and  cloud 
base  was  at  1.2  km.  A  cloud  which  developed  over  the 
NASA  field  mill  network  was  subsequently  penetrated  by 
the  T-28  at  an  altitude  of  4.8  km  msl.  The  CP-2  radar 
captured  the  full  cloud  evolution  from  first  echo  through 
the  mature  stage.  The  cloud  was  first  detected  (8.5  dBZ 
at  4.3  km  height,  range— 23.5  km)  at  18:21  UTC  (hence¬ 
forth  all  times  are  UTC  and  all  heights  are  msl).  For 
a  few  minutes  the  dual-frequency  (S/X-band)  reflectivity 
ratio,  DFR,  was  8  -  9  dB  giving  evidence  of  partial  Bragg 
scattering,  the  Zdr  being  0  dB.  At  18:25,  Zh  increased  to 
10  dBZ,  the  DFR  decreased  to  0  dB  while  Zdr  jumped 
to  1.7  dB  at  5  km  height.  This,  we  believe,  is  strong  ev¬ 
idence  for  raindrop  formation  by  the  warm  rain  process. 
At  18:28  a  few  big  drops  (Zdr  ~  3  -  4  dB)  were  detected 
falling  to  the  ground  while  an  anvil  from  a  mature  storm 
(about  10  km  to  the  SE)  extended  over  the  cloud  of  in¬ 
terest  at  9  km  height.  It  appears  from  the  radar  data 
that  ice  crystals  from  the  anvil  descended  into  the  upper 
part  of  this  cloud,  freezing  it.  Zdr  was  nearly  0  dB  above 
5.5  km,  with  large  drops  continuing  to  fall  out  from  the 
cloud.  The  reflectivity  aloft  (6  km)  intensified  to  48  dBZ 
at  18:36  with  Zdr  being  0  dB,  indicating  graupel  forma¬ 
tion.  We  refer  to  this  cloud  as  cell  A,  to  distinguish  it 
from  new  growth  that  occurred  on  the  NW  side  of  cell 
A,  a  region  with  persistent  updraft.  The  first  evidence  of 
new  growth  was  at  18:35  and  this  cloud  will  be  called  cell 
B  as  its  evolution  could  be  separately  tracked.  Its  micro¬ 
physical  evolution  was  rapid  and  very  different  from  cell 
A  and  was  responsible,  we  believe,  for  cloud  electrification 
which  started  at  18:41  with  the  first  CG  strike  at  18:43. 

2.  CLOUD  EVOLUTION 

Fig.l  shows  CAPPIs  at  5.5  km  height  with  T-28  tracks 
which  penetrated  the  cloud  at  18:36  and  18:41  overlaid. 
The  grey  scales  depict  Zdr,  the  darker  shades  represent¬ 
ing  larger  values.  Fig. la  shows  cell  A  about  12  min  after 
first  echo.  Fig.lb  shows  a  positive  Zdr  column  (cell  B) 
while  cell  A  has  glaciated  (~  0  dB  Zdr).  Fig. 2  shows  ver¬ 
tical  sections  taken  along  the  shear  direction  (SE-NW); 


Fig.l  CAPPI  at  5.5  km  height  of  Zh  contours  (start  from  10 
dBZ,  step  10  dB)  with  Zir  grey  scales  (the  intervals 
being  0.5,1.5,3,4  dB). 
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Fig. 2  Vertical  sections  along  an  SE-NW  line  in  Fig.l.  X=0 
here  is  the  center  of  the  reflectivity  core  in  Fig.l.  Zdr 
grey  scales  same  as  in  Fig.l. 
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positive  (negative)  X  values  represent  locations  SE  (NW) 
of  the  main  echo.  Note  the  positive  Zdr  column  in  Fig.2b 
(cell  B).  The  height  profile  of  Zdr  represents  updraft  con¬ 
ditions  in  cell  B  and  downdraft  in  cell  A.  Fig.  3  sum¬ 
marizes  T-28  data  for  the  18:36  (solid  lines)  and  18:41 
(dashed  lines)  penetrations,  which  were  respectively  along 
the  SE  — *  NW  and  NW  — •  SE  directions.  At  18:36  the 
T-28  detected  a  peak  downdraft  of  9  m/s  just  below  the 
main  echo  in  Fig.2a  (cell  A).  The  PMS  2D-P  probe  as 
well  as  the  Cannon  camera  showed  a  mixture  of  very  wet 
graupel  and  drops  in  the  downdraft  with  maximum  sizes 
of  6  -  7  mm  as  shown  in  Fig. 3.  The  LDR  was  about  - 
23  dB  in  the  downdraft,  consistent  with  wet  graupel.  A 
peak  updraft  of  13  m/s  (LWC  ~  1.5  g/m3)  was  detected 
about  2.5  km  NW  of  the  main  echo  in  Fig.la.  This  was 
well  correlated  with  the  first  radar  detection  of  cell  B  at 
18:36  at  5  km  height.  The  T-28  pilot  reported  moder¬ 
ate  rain  early  during  penetration  of  the  updraft.  Neither 
the  P-probe  nor  the  Cannon  camera  detected  any  drops 
in  the  updraft;  however,  the  foil  impact  or  registered  a 
few  circular  impressions  (drops)  of  1-2  mm  diameter  in 
the  updraft.  On  the  return  penetration  at  18:41,  the  T- 
28  sampled  cell  B  (updraft  peak  ~  15  m/s,  LWC  ~  2 
g/m3)  and  cell  A  (downdraft  peak  ~  7  m/s,  wet  grau¬ 
pel  6  -  7  mm).  Within  the  updraft  (positive  Zdr  column, 
Fig.2b),  the  P-probe  detected  a  few  large  drops  (5-6 
mm)  while  the  foil  impactor  registered  eight  circular  im¬ 
pressions,  the  largest  being  2.5  mm,  and  about  10s  later 
a  large  drop  (~  6  mm).  The  T-28  pilot  reported  moder¬ 
ate  rain  in  the  updraft.  These  data  are  consistent  with 
the  positive  Zdr  column  in  Fig.2b  (cell  B)  consisting  of  a 
low  concentration  of  large  drops.  The  pilot  also  reported 
lightning  during  this  penetration  as  well  as  a  transition 
from  rain  to  ice  on  the  windscreen  as  he  left  the  updraft. 
A  dry /wet  growth  modeling  calculation  showed  that  the 
larger  graupel  (6-7  mm)  could  be  in  wet  growth.  Peak 
LDR  and  X-band  specific  attenuation  values  of  -18  dB 
and  1  dB/km  were  observed  at  the  T-28  flight  altitude  in 
the  downdraft,  which  are  also  consistent  with  wet  graupel. 

Fig. 4a  shows  the  time-height  reflectivity  profile  of  cell 
A  starting  at  18:23;  contours  start  at  0  dBZ  and  increment 
by  5  dB-  The  dark  solid  line  shows  the  maximum  height 
reached  by  the  2  dB  Zdr  contour.  Note  the  overlying  anvil 
from  a  mature  storm  about  10  km  to  the  5E  referred  to 
earlier.  The  arrows  represent  peak  vertical  winds  from  the 
T-28  penetrations  at  18:36,  18:41  and  18:48  and  one  pen¬ 
etration  of  the  NOAA  P-3  at  18:45.  The  top  panel  shows 
the  vertical  electric  held  from  a  surface  field  mill  located 
under  the  cloud  at  X=13,  Y=-24  km  (see  Fig.lb).  The 
dark  vertical  bars  are  the  maximum  vertical  component 
of  the  E-field  from  the  T-28  field  mills.  Horizontal  dark 
bars  represent  the  duration  of  the  penetration.  The  solid 
triangle  at  1200s  is  the  time  of  the  first  cloud- to-ground 
lightning  detected  by  the  NASA  LLP  network  (18:43:27  at 
X=13.5,  Y=-24.6  km).  Before  600s  the  surface  field  mill 
activity  is  due  to  a  strong  storm  located  10  km  W/SW. 
Fields  with  fair  weather  polarity  start  at  700s  with  a  re¬ 
versal  evident  at  900s.  At  1050s  (18:41  T-28  penetration 
in  Fig.3)  the  electric  fields  build  up  at  the  T-28  altitude 
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Fig.3  Data  from  T-28  showing  up/downdraft  (velocity),  liq¬ 
uid  water  content  (lAV c),  vertical  E-iield  and  maximum 
partical  size.  The  time  t=0  is  the  starting  time  for  each 
penetration. 


with  the  CG  strike  following  at  1200s.  As  noted,  the  T- 
28  pilot  reported  lightning  during  the  18:41  penetration, 
the  surface  field  mill  shows  numerous  strikes  after  this 
time.  Cell  A  shows  a  general  sinking  of  the  reflectivity 
with  time  after  700s  (~  18:34)  with  some  indication  of 
pulsation.  The  maximum  height  of  the  2  dB  Zdr  con¬ 
tour  also  sinks  slowly,  correlating  well  with  the  persistent 
downdraft  measured  by  the  two  aircraft. 

The  evolution  of  cell  B  in  Fig.4b  is  quite  different.  It 
was  first  detected  as  new  growth  at  7G0s.  The  strong 
updraft  detected  by  the  T-28  (Fig.3)  correlates  well  with 
the  positive  Zdr  column,  the  2  dB  Zdr  height  rising  to 
6  km  (— 8*C).  Subsequently,  following  what  appears  to 
be  freezing  of  the  raindrops  the  Zdr  height  falls  while  the 
cloud  top  grows  substantially,  possibly  due  to  the  associ¬ 
ated  latent  heat  release.  By  the  time  of  the  T-28’s  18:41 
penetration,  the  50  dBZ  core  is  at  6.5  km.  This  was  also 
a  region  of  wet  growth  as  predicted  by  the  growth  model. 
The  2  dB  Zdr  height  doesn’t  decrease,  remaining  near 
the  —  3*C  level  indicating  persistent  updraft  as  corrobo¬ 
rated  by  the  P-3  and  T-28  penetrations  at  18:45  (1350s) 
and  18:48  (1500s),  respectively.  The  surface  field  mill  lo¬ 
cated  near  cell  B  at  X=7,  Y =-22  km  (see  Fig.l)  shows 
fair  weather  polarity  until  1050s  (~  18:40),  with  subse¬ 
quent  electrification  until  the  CG  strikes  are  indicated 
(top  panel).  The  T-28  field  mills  also  show  electrification 
at  the  time.  It  appears  that  mixed  phase  microphysics 
plus  the  up/downdraft  structure  played  important  roles 
in  electrifying  the  clouds.  Detailed  analysis  of  the  aircraft 
data  and  other  radar  signatures  such  as  LDR  and  Aj  will 
add  important  details  to  the  cloud  evolution. 
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Fig.4  Time-height  contours  of  for  cell  A  (a)  end  cell  B  (b).  Dark  solid  line  is  the  maximum  height  of  the  2  dB  Z4, 
level  reached  at  various  times.  Vertical  arrows  are  peak  updraft  speeds  from  T-28  (3)  and  P-3  (1).  Top  panels  show 
maximum  vertical  El-field  from  T-28,  £«  (dark  solid  bars),  and  the  El-field  from  surface  field  mill  located  near  the 
cells,  E,  (light  continous  line).  The  solid  triangle  mark  in  (a)  at  1210  s,  and  (b)  at  1260  s,  is  the  first  CG  strike  from 
the  cell. 
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1.  INTRODUCTION 

The  CaPE  project  offered  the  first  opportunity  to  col¬ 
lect  dual-polarized,  dual-frequency  radar  data  in  the  sub¬ 
tropical  regime  of  centra^  Florida.  This  paper  is  based 
on  CP-2  time  series  data  collected  at  vertical  incidence 
in  light  rain,  and  in  a  squall  line  case  at  close  range. 
Each  time  series  record  consisted  of  128  H  and  128  V 
samples  at  S-band,  and  256  H  samples  at  X-band.  From 
these  data,  reflectivity  at  horizontal  polarization  2a,  Z&r, 
copolar  correlation  coefficient  pa.(O),  differential  propa¬ 
gation  phase  <bdP  and  dual-frequency  (S/X)  reflectivity  ra¬ 
tio  DFR  were  calculated.  A  new  algorithm  for  calculating 
specific  attenuation  As  at  X-band  and  specific  differential 
phase  Kip  is  used  based  on  Hubbert  et  al.  (1993),  which 
accounts  for  hail  signal  and  phase  shift  upon  backscat- 
ter  fluctuations  which  may  be  superimposed  on  the  range 
profiles  of  DFR  and  4>dp,  respectively. 

2.  CP-2  RADAR  MEASUREMENTS 

The  radar  data  reported  here  were  collected  on  24  Au¬ 
gust,  1991  of  the  CaPE  project. 

2.1.  Vertically  Pointing  Data 

Raindrops  at  vertical  incidence  form  a  standard  target 
for  polarization  measurements  because  of  rotational  sym¬ 
metry  of  the  drops  when  viewed  vertically.  In  addition,  if 
the  antenna  is  rotated  a  full  360*,  the  polarization  plane 
also  rotates  with  it.  For  this  method  to  work  the  melting 
level  must  be  located  beyond  the  far-field  of  the  antenna 
(in  this  case  ~  1.5  km),  a  common  situation  in  Florida 
but  not  in  Colorado,  for  example. 

Fig.l  shows  vertical  profiles  of  Z/,,  phv(0),<t>dP,  and 
DFR  in  light  rain  at  90*  elevation  angle.  These  data 
have  not  been  filtered  in  range.  The  vertical  bars  in  each 
case  show  the  sample  standard  deviation.  The  Zs  bright- 
band  is  clearly  visible  at  3.5  km  (all  heights  are  agl),The 
dip  in  p*.(0)  to  0.92  at  3.25  km  is  also  clearly  noted. 
Note  that  pa.( 0)  is  estimated  from  the  power  time  se¬ 
ries  (see  Liu  et  al.,  these  Proceedings).  The  mean  value 
of  psv(0)  in  light  rain  below  the  bright-band  is  around 
0.982.  It  is  well  known  that  pa«(0)  should  be  very  close 
to  unity  in  light  rain  even  at  horizontal  incidence  with 
values  of  0.997  quoted  by  Illingworth  &  Caylor  (1991). 
At  vertical  incidence  the  deviation  from  1  should  be  even 
less.  One  possible  reason  for  the  measured  ph»(0)  being 
so  low  is  that  the  phase  difference  between  the  HH  and 
VV  antenna  patterns  within  the  main  lobe. deviates  due 
to  feedhorn  imperfections  or  loss  of  precision  of  the  reflec¬ 
tor  shape.  Simulations  show  that  this  effect  is  sufficient 
to  drop  p/i„(0)  to  0.99  (see  Xiao  et  al.,  these  Proceed- 
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Fig.]  Range  profiles  of  mean  and  standard  deviation  of 

(a)  and  Pa.(O), 

(b)  <t>dp  and  DFR  front  vertical  pointing  data. 


ings).  The  increase  in  a  with  height  >  5.5  km  is  due  to 
SNR  becoming  smaller.  The  large  a  at  the  p*„(0)  dip 
is  consistent  with  theory  (see  Liu  et  al.,  these  Proceed¬ 
ings)  and  with  the  results  of  Illingworth  &  Caylor  (1991). 
The  system  <bdp(~  —50.5®)  can  be  obtained  from  Fig. lb 
since  at  vertical  incidence  no  differential  propagation  ef¬ 
fects  are  expected  in  light  rain.  The  DFR  results  are  also 
shown  in  Fig.lb.  Note  the  positive  fluctuation  of  about 
1.5  dB  in  the  mean  at  3.75  km,  very  close  to  the  peak 
Zh.  This  is  suggestive  of  Mie  scattering  effects  due  to 
large,  wet  snowflakes  with  maximum  size  around  10  min. 
An  assessment  of  the  overall  accuracy  in  derived  specific 
attenuation  (Aj.dBkm"1)  and  specific  differential  phase 
(Kip.'  km~l)  can  be  made  since  the  mean  should  equal 
zero  in  light  rain.  Errors  may  be  due  to  statistical  fluctua¬ 
tions,  system  effects,  and  the  algorithms  used  to  estimate 
As  and  KdP-  Here  we  use  an  adaptation  of  range  filters 
described  by  Hubberltt  al  (1993)  to  calculate  A3  and  AuP. 
The  accuracy  in  estimating  the  mean  value  is  defined  here 
as  the  sum  of  the  absolute  value  of  the  bias  and  the  95% 
confidence  interval  which  is  0.2  -  0.25  (dBkm~l or* km~l ) 
for  both  A3  and  Kip.  Note  the  condition  at  vertical  in¬ 
cidence,  i.e.,  high  SNR,  narrow  spectrum  width  and  very 

high  Ph.(0). 

2.2.  Rainfall  Data 
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Fig.2  (a)  Zh  contours  with  Zdr  grey  scales  overl&yed; 

(b)  mean  and  95%  C.I.  of  ph»( 0)  vs  Zdr', 

(c)  mean  and  95%  C.I.  of  Aj  vs  Kdpi 

(d)  mean  and  95%  c.l.  of  Kdp  vs  Z*. 


On  24  August  1991  &  squall-tine  approached  the  radar 
from  the  W-NW.  Zh,  Zdr,  DFR,  Qdp  and  have 

been  filtered  in  range  using  an  HR  filter  which  sharply 
attenuates  spatial  fluctuations  250  m  and  less. 

Fig. 2a  shows  a  PP1  scan  of  Zh  at  1.5*  elevation  an¬ 
gle  derived  from  the  central  Florida  data.  Grey  scales 
correspond  to  Zdr  with  darker  shades  representing  higher 
values. 

Fig.2b  shows  a  plot  of  phv(O)  vs  Zjr  obtained  from  the 
1.5*  elevation  PP1  in  Fig. 2a  between  38  -  48  km.  The 
mean  p*„( 0)  does  decrease  with  increasing  Zdr  in  agree¬ 
ment  with  theory,  but  the  p*,v(0)  values  are  significantly 
lower  compared  to  both  theory  and  the  experimental  re¬ 
sults  of  Illingworth  &  Caylor  (1991).  We  refer  to  A'iao  et 
al.  (these  Proceedings)  for  a  further  discussion  of  antenna 
effects  related  to  the  lowered  values  of  pa«(0). 

Fig.2c  shows  Aj  vs  Kdp  in  rain  from  data  collected  at 
1.5*  elevation  (37  -  48  km)  and  subsequently,  at  elevation 
angles  2  -  4*  (14  -  26  km).  The  mean  and  95%  confidence 
intervals  are  shown  in  Fig. 2c  and  are  marked  as  data  from 
central  Florida.  Another  curve  marked  ‘Colorado’  is  CP- 
2  time  series  data  from  a  rainshaft  in  a  Colorado  storm 
at  low  elevation  angle  (1*)  at  range  82  -  93  km.  These 
data  have  been  processed  exactly  the  same  as  for  the  cen¬ 
tral  Florida  data.  In  the  Kdp  range  1  —  2.25*krn_,)  the 
two  data  sets  have  different  slopes;  the  central  Florida 
data  has  a  slope  of  1.05  while  the  Colorado  data  has 
a  slope  of  0.78.  Also,  for  a  given  Kdp,  the  Aj  values 
in  Colorado  are  consistently  less  than  in  central  Florida. 
We  believe  that  these  features  are  statistically  significant 
and  are  a  manifestation  of  microphysical  differences  be¬ 
tween  the  two  regimes,  principally  ralated  to  drop  size 
distributions,  and  secondarily  to  mean  raindrop  shape  vs 
size  relations.  The  Colorado  slope  of  0.78  is  very  close 
to  that  predicted  theoretically  by  Bringi  et  al.  (1990) 
for  Gamma  DSDs  (0.8),  and  for  experimental  DSDs  (col¬ 
lected  in  Norman,  OK)  from  surface  distroineter  (0.79) 
assuming  equilibrium  raindrops  shape  vs  size  relation  due 
to  Green  (1975).  The  results  from  Fig. 2c  are  independent 
of  absolute  system  gain  and  as  a  consequence  are  more 
reliable. 

Fig. 2d  shows  the  mean  Kdp  vs  Zh  for  the  central 
Florida  and  Colorado  datasets.  The  solid  line  is  the 
mean  relationship  from  simulations  using  Gamma  DSDs. 
Again,  for  a  given  Kdp,  the  Zh  in  Colorado  is  as  expected 
larger  than  in  central  Floridais. 

2.3.  Vertical  Cross  Sections  through  A  Squall  Line 

We  now  show  vertical  cross  sections  through  the  most 
intense  part  of  the  24  August  squall  line.  Fig.3  shows  an 
RHI  scan  taken  to  the  west  of  the  radar.  Darker  shades 
represent  higher  values  of  Zdr  and  As.  Contours  start  at 
10  dBz  and  increament  by  10  dBz.  These  data  were  col¬ 
lected  in  the  normal  operating  mode  of  the  CP-2  radar 
which  does  not  allow  for  computation  of  differential  phase 
and  ph »(0).  Fig. 3a  shows  Zh  contours  with  Zdr  grey  scales 
overlayed,  while  Fig.3b  show  contours  of  Zh  with  Aa  over- 
layed.  The  leading  edge  of  the  squall  line  is  at  18  km 
range.  Note  the  positive  Zdr  column  at  22  km,  with  the 


Fig.3  Zt,  contours  with  (a)  Zdr,  (b)  A3 .  grey  scales  overlayed. 


2  dB  contour  reaching  6.5  km  agl.  The  peak  Zdr  is  4.5 
dB  at  5  km  agl.  The  A3  core  is  located  at  4.5  km  agl  in 
Fig. 3b  in  the  range  sector  24  -  27  km  where  Zdr  is  1.5 
-  2.0  dB.  This  >13  core  represents  the  main  precipitation 
mass.  On  the  inflow  side,  a  secondary  A3  maximum  can 
be  located  at  6  km  agl  in  the  range  22  -  23  km  which  also 
coincides  with  high  Zh  (60  dBZ)  and  high  Zdr  (3.5  -  4 
dB).  This  region  consists  of  large  raindrops  suspended  by 
the  updraft. 

About  3  minutes  later,  time  series  data  were  collected 
through  the  same  storm  at  the  same  location  as  described 
in  Fig.3.  Range  profiles  at  elevation  angle  11’  are  shown 
in  Fig. 4.  This  profile  intersects  the  main  precipitation 
mass  shown  in  Fig.3b  at  4.7  km  agl,  range  =  25  km. 
Fig.4a,b  shows  range  profiles  of  Zh,  Zdr,  4>dp  and  Kdp. 
Note  the  peak  Kdp  of  4.75*Jtm-1  corresponding  to  a  rain¬ 
fall  rate  of  150mmA~‘ .  Note  the  decrease  in  Zdr  from  3.75 
dB  at  23.5  km  to  0.5  dB  at  29  km,  indicating  the  transi¬ 
tion  from  rain  to  ice.  Fig.4c,d  show  range  profiles  of  Zh, 
Phv( 0),  DFR  and  A3.  From  23  to  29  km,  Ph*(0)  (unfil¬ 
tered  in  range)  decreases  from  0.98  to  0.93  indicating  a 
transition  from  rain  to  ice.  Beyond  30  km,  the  average 
Phv(0)  is  0.96,  and  beyond  35  km  it  is  0.985.  The  peak  A3 
is  4.7  dBkm~l  and  its  location  coincides  with  the  peak 
Kdp  location.  The  DFR  curve  decreases  with  range  from 
30  km  onwards  as  the  X-band  signal  falls  below  noise  level 
due  to  attenuation. 

3.  CONCLUSIONS 

Vertical  incidence  data  in  light  rain  revealed  a  dip  in 


Fig.4  Range  profiles  of  (a)  Zh  and  Zrfr;  (b)  4>dp  and  Kdp\ 
(c)  Zh  and  pi ,,(0);  (d)DF/t  and  A3. 


pi ,*(0)  just  below  the  bright-band,  as  well  as  weak  Mie 
scattering  effects.  The  accuracy  in  estimating  Aj  and 
Kdp  (from  the  vertical  incidence  data)  was  established  to 
be  0.2  dBkm~l  or  mkm~l. 

Data  from  central  Florida  and  Colorado  were  com¬ 
pared  in  convective  rainfall  which  showed  differences  at¬ 
tributable  to  drop  size  distribution  variations  and  to  a 
lesser  extent  variations  in  the  mean  raindrop  shape  vs 
size  relation. 

Vertical  sections  of  data  through  the  intense  portion  of 
a  squall  line  revealed  interesting  tnicrophysical  features 
such  as  a  positive  Zdr  column  adjacent  to  the  leading 
edge,  the  transition  from  rain  to  ‘rain  +  ice’  to  ice  as 
revealed  by  psv(0),  and  the  large  values  of  A3  and  Kdp 
within  the  main  precipitation  mass. 
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MICROPHYSICAL  AND  ELECTRICAL  EVOLUTION  OF  A  CONVECTIVE  STORM  USING 
MULTIPARAMETER  RADAR  AND  AIRCRAFT  DATA  DURING  CaPE 
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1.  INTRODUCTION 

On  9  August  1991  of  the  CaPE  project  the  transition 
from  early  echo  to  mature  stage  was  captured  by  the  CP-2 
radar  as  well  as  two  King  Airs  (NCAR  KA  and  Wyoming 
KA),  the  NOAA  P-3  and  the  T-28  aircrafts.  It  was  the 
only  such  case  in  CaPE  involving  good  coordination  be¬ 
tween  the  multiple  aircraft  and  radar  permitting  a  study 
of  the  microphysical  and  electrical  evolution  of  the  par¬ 
ticular  cloud.  A  companion  paper  by  Bringie t  al.  (these 
Proceedings)  deals  with  the  evolution  of  another  adjacent 
cloud  (about  7  km  away)  which  was  penetrated  by  the 
T-28.  The  general  environmental  conditions  were  quite 
similar  for  both  clouds  and  we  refer  to  that  paper  for 
details.  However,  the  microphysical  evolution  has  some 
differences  between  the  two  clouds  that  are  interesting. 

The  cloud  was  first  detected  (peak  of  3  dBZ  at  height 
=  2.6  km  and  range  =  25  km)  at  17:48  UTC  (hence¬ 
forth  all  height  are  msl  and  time  are  UTC).  The  dual- 
frtquency  (S/X-band)  reflectivity  ratio,  DFR,  was  15  dB 
indicating  dominance  by  Bragg  scattering.  At  17:51,  the 
peak  reflectivity  increased  to  22.5  dBZ,  DFR  decreased 
to  0  dB  and  Zdr  increased  to  0.75  dB,  evidence  of  co¬ 
alescence  formation  of  drops.  About  10  minutes  later 
(18:00)  a  strong  positive  Zdr  column  had  formed,  simi¬ 
lar  to  those  reported  by  Illingworth  (1992)  in  Alabama 
during  MIST.  The  warm  cloud  with  a  positive  Zdr  col¬ 
umn  persisted  until  18:09  and  then  the  cloud  glaciated, 
releasing  latent  heat  and  accompanied  by  strong  vertical 
growth.  A  weak  E-field  (~  5  kV/m)  detected  by  NCAR 
KA  at  18:12  forced  it  away  from  the  cloud.  An  enhanced 
E-fleld  was  detected  by  the  NOAA  P-3  at  18:16.  The  first 
cloud-to-ground  lightning  from  the  cloud  was  detected  by 
the  NASA  LLP  network  at  18:18. 

2.  CLOUD  EVOLUTION 

Surface  winds  were  from  the  W/NW  near  the  Cape. 
The  15:00  sounding  at  Cape  Canaveral  indicated  light 
westerly  winds  up  to  6  km  turning  E/SE  above  7.5  km. 
The  0*C  level  was  at  4.6  km,  with  cloud  base  near  1.2 
km.  The  LWC  observations  from  aircraft  were  less  than  1 
g/m3  indicating  a  large  amount  of  mixing  in  these  clouds 
(adiabatic  LWC  from  sounding  was  ~  6  g/m3). 

The  cloud  (A)  of  interest  was  detected  (17:55)  about  4 
km  NW  of  another  cloud  (B)  that  was  in  a  more  mature 
phase.  Both  clouds  were  about  10  km  NE  of  an  intense 
convective  line  of  thunderstorms.  At  18:01,  cloud  B’s  re¬ 
flectivity  core  (~  47  dBZ  at  4.5  km)  started  to  descend 
while  the  now  glaciated  upper  part  of  the  cloud  (4-7  km) 
began  to  tilt  to  the  NW  in  response  to  E/SE  winds  aloft. 


Fig.l  CAPPI  at  5.5  km  height  of  Zh  contours  with  Zdr  grey 
scales.  Darker  shades  are  larger  ZdT,  the  intervals  being 
.5, 1.5,3, 4  dB. 


At  18:04,  the  warm  cloud  A  which  was  on  the  up-shear 
side  (top  ~  5  km)  had  the  tilted  upper  part  of  cloud  B 
overhead.  It  is  possible  that  surface  outflow  from  cloud 
B  enhanced  low  level  convergence  near  cloud  A.  It  is  also 
possible  that  tiny  ice  crystals  or  graupel  from  cloud  B 
may  have  descended  and  melted  into  cloud  A  forming 
tiny  drops  that  can  grow  quickly  to  very  large  raindrops. 
Indeed  at  18:08  the  3  dB  Zdr  contour  in  cloud  A  reached 
up  to  6  km,  with  maximum  Zdr  of  4  dB  at  4  km  height 
(prior  to  this  the  maximum  Zdr  was  2.5  dB  at  3  km,  and 
cloud  top  was  5  km).  Also,  cloud  B  loses  its  identity  after 
18:05,  more  or  less  merging  with  cloud  A. 

Fig. la, b  show  CAPPI  of  cloud  A  at  5.5  km  height  at 
(a)  18:08  and  (b)  18:15;  contours  of  Zh  start  at  10  dBZ 
and  increment  by  10  dB  while  Zdr  is  shown  as  a  grey  scale 
with  darker  shades  representing  larger  Zdr-  Cloud  A  is 
centered  at  X=13,  Y=-27  km  in  Fig.la  which  also  shows 
the  positive  Zdr  column.  In  Fig. lb  the  strength  of  the 
column  has  weakened  considerably.  The  NOAA  P-3  pene¬ 
trations  at  18:09  and  18:15  are  shown  in  Fig.l.  Figs.2a,b,c 
show  vertical  sections  taken  along  a  line  oriented  SE-NW 
(132*  from  N)  in  Fig.l;  positive  (negative)  X  represents 
distance  SE  (NW)  from  the  core.  The  positive  Zdr  col¬ 
umn  in  Fig.2a  is  clearly  visible  at  18:08,  while  at  18:12 
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Fig. 2  Vertical  sections  along  an  SE  -  NW  line  in  Fig.l.  X=0 
here  is  the  center  of  the  core  in  Fig.l.  Zdr  grey  scales 
same  as  in  Fig.  1 . 


(Fig. 2b)  the  cloud  has  glaciated.  At  18:15  the  reflectivity 
core  was  descending  and  the  Zdr  structure  conforms  to 
the  more  usual  ice-melting-drops  situation. 

The  NOAA  P-3  made  its  penetrations  at  am  altitude  of 
6.4  1cm  (— 10*C).  Fig.3  shows  some  of  the  aircraft  mea¬ 
sured  data  during  the  18:09  (solid  line)  and  18:15  (dashed 
line)  penetrations.  The  top  panel  shows  up/downdraft 
and  liquid  water  content  (LWC),  while  the  bottom  panel 
shows  the  vertical  E-field.  At  18:09  the  P-3  (going 
SE  — *  NW)  was  flying  on  top  of  the  positive  Zdr  column 
(shown  in  Fig.2a).  The  main  precipitation  core  (Zt>  ~  50 
dBZ)  was  located  about  1  km  below  the  flight  altitude. 
Peak  Zdr  (~  3  dB)  was  located  at  4.4  km.  The  updraft 
location  coincides  well  with  the  positive  Zdr  column.  Be¬ 
fore  encountering  the  updraft  (from  the  SE  side  of  core) 
the  2D-C  probe  showed  rough  edged  graupel  (maximum 
dimension  ~  2  mm).  Upon  entering  the  updraft,  the  im¬ 
ages  were  very  smooth  with  particle  types  being  drops 
and  wet  graupel.  A  3  mm  raindrop  was  detected  in  the 
center  of  the  second  updraft  peak  (~  70s  from  0  in  Fig.3). 
Within  the  updraft  region  a  high  LDR  layer  was  detected 


with  values  around  -17  dB  with  X-band  specific  attenua¬ 
tion  (Aj)  of  1  dB/km.  Zdr  was  around  1-1.5  dB.  The 
large  LDR  signature  indicates  wet  graupel  while  both  Zdr 
and  Aj  are  indicative  of  supercooled  drops  consistent  with 
the  2D-C  images.  The  E-field  is  very  weak  at  18:09. 

The  NCAR  and  Wyoming  King  Air  (NKA,  WKA)  each 
made  a  NW  —  SE  penetration  (-S’C^.S’C)  centered 
at  18:11.  At  this  time  (Fig.2b)  the  cloud  glaciated  with 
cloud  top  rising  from  7.5  km  (at  18:08)  to  9.0  km  (at 
18:12).  The  NKA’s  penetration  at  5.5  km  went  just  on 
top  of  the  column  while  the  WKA  was  in  the  high  Zdr 
region  (3  -  3.5  dB)  at  4  km  height  (Fig. 2b).  The  NKA 
detected  a  peak  updraft  of  8  m/s  on  top  of  the  column 
and  downdraft  peak  of  -3  m/s  adjacent  to  it.  At  this 
flight  level,  a  layer  of  high  LDR  (-18  to  -20  dB)  and  A3 
(1  -  1.5  dB/km)  was  noted  at  the  up/down  transition. 
The  NKA  detected  fields  of  5  -  6  kV/m  and  soon  exited 
the  cloud.  The  WKA  at  lower  altitude  (4  km)  detected 
a  peak  updraft  of  9  m/s  followed  by  a  peak  downdraft  of 
-4  m/s  separated  by  900  m.  Large  drops  (5-7  mm)  were 
detected  by  the  P-probe  during  this  penetration. 

The  P-3  made  its  reverse  penetration  (going  NW  — » 
SE)  at  18:15  with  data  shown  in  Fig.3.  The  vertical  ve¬ 
locity  has  decreased  substantially  and  the  El- fields  have  in¬ 
creased  dramatically.  Fig.2c  shows  the  vertical  Zdr  struc¬ 
ture.  At  the  P-3  flight  level  (6.4  km),  Zdr  is  nearly  0  dB 
while  LDR  is  -26  dB  and  A3  is  near  zero.  The  storm  is 
in  its  mature  phase  and  the  main  precipitation  core  (Zi, 
peak  of  57  dBZ)  is  descending  to  the  surface.  The  2D-C 
probe  detects  a  much  higher  concentration  of  rough-edged 
graupel  (maximum  sizes  ~  2  -  3  mm). 

Fig. 4  shows  a  time-height  cross  section  of  cloud  evo¬ 
lution  from  first  echo  (17:48:00  is  t=0)  through  a  30 
min  period.  The  dark  solid  line  is  the  maximum  height 
reached  by  the  2  dB  Zdr  contour.  Note  the  maximum 
height  of  6  km  reached  at  ~  18:08  (1250s  from  t=0)  fol¬ 
lowing  which  there  is  a  steady  decrease  to  2  km.  The 
vertical  arrows  represent  updraft  detected  by  the  P-3 
(6.5  km),  NKA  (5.5  km)  and  T-28  (4.8  km)  at  various 
times  in  the  cloud  evolution.  The  updraft  magnitude  de¬ 
creases  with  time  are  roughly  correlated  with  the  maxi¬ 
mum  height  of  the  2  dB  Zdr  contour.  The  top  panel  shows 
the  detection  of  the  first  cloud-to-ground  lightning  by  the 
NASA  LLP  network  at  18:18:12  at  X=14.1,  Y=-29.5  km 
(see  Fig. lb).  It  appears  that  mixed  phase  microphysics 
plus  the  up/downdraft  structure  played  important  roles 
in  electrifying  the  cloud.  Detailed  forthcoming  analysis 
of  all  the  aircraft  data  with  radar  LDR  and  A3  signa¬ 
tures  will  add  important  details  to  cloud  microphysical 
and  electrical  evolution.  This  study  can  be  compared  to 
another  cloud  evolution  described  by  Bringi  et  al.  (these 
Proceedings). 
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Data  from  NO  A  A  P-3  showing  up /downdraft  (velocity),  liquid  water  content  (LWC)  and  vertical  E-field.  The 
t=0  is  the  starting  time  for  each  penetration. 
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Fig.4  Time-height  contour*  of  Z) Dark  solid  line  is  the  maximum  height  of  the  2  dB  Z^r  level  reached  at  various  times. 
Vertical  arrows  are  updraft  speeds  from  various  aircraft.  Top  panel  shows  vertical  El-field  from  P-3.  The  solid  triangle 
mark  at  1840s  is  the  first  CG  strike. 
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